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ABSTRACT 


Problems  associated  with  the  realization  of  a  high  energy, 
highly  coherent  Raman  laser  have  been  investigated.  Since  the  dis¬ 
covery  that  the  self-focusing  of  pump  radiation  and  scattered  radi¬ 
ation  is  present  in  most  laboratory  Raman  lasers  and  tends  to  spoil 
the  coherence  of  the  output,  the  nature  and  mechanisms  of  this  self¬ 
trapping  have  been  the  primary  objects  of  study  in  order  that  it 
might  be  prevented.  Time  resolved  photographic  studies  have  been 
made  which  show  that  the  small  filaments  of  self- trapped  light  in 
the  pump  beam  persist  for  only  a  small  fraction  of  the  beam  pulse 
duration  and  are  sufficiently  intense  to  resolve  the  apparent  dis¬ 
crepancy  between  the  average  intensity  and  observed  Raman  gains. 
Time  integrated  photographic  studies  and  photometric  studies  showed 
that  only  purely  linearly  polarized  light  preserves  its  polarization 
after  self-trapping  in  CS£  and  nitrobenzene.  Attempts  to  compare 
measured  nonlinear  indices  for  both  circular  and  linear  polari¬ 
zations  with  those  calculated  for  various  mechanisms  were  frustrated 
by  this  instability  in  the  propagation  of  circularly  polarized  light. 

The  instability  was  explained  qualitatively  and  indicates  that  electro  - 
striction  is  not  a  dominant  mechanism.  A  classical  statistical 
mechanical  calculation  of  the  nonlinear  index  of  a  liquid  (with  no 
macroscopic  density  changes)  showed  that  the  local  redistribution 
of  molecules  is  generally  important.  This  effect  will  cause  perfectly 
symmetric  molecules  to  have  a  significant  nonlinear  index.  It  is 
therefore  concluded  that  it  is  very  unlikely  that  any  room  temperature 
liquid  will  perform  as  a  high  energy  Raman  laser  material  without 
beam  deterioration  due  to  self-focusing. 


v 


I. 


INTRODUCTION 


The  growth  of  our  understanding  of  stimulated  Raman  scattering 
(SRS)  can  be  characterized  as  a  series  of  discoveries,  preceded  and 
followed  by  periods  of  wide  divergence  of  theory  and  experiments.  It 
could  be  said  that  our  progress  has  been  one  of  discovering  one  prob¬ 
lem  after  another  and  slowly  resolving  these  problems  by  discovery  of 
new  effects.  We  list  here  some  of  the  problems  which  have  plagued 
the  development  of  SRS;  some  of  them  have  not  been  resolved  within 
the  period  covered  by  this  report. 

The  first  problem  to  attract  attention  concerned  the  incorrect 
emission  angles  of  the  anti-Stokes  lines;  the  second  problem  involved 
the  large  spectral  broadening  due  to  multimode  pumping,  and  the  third 
concerned  the  large  amount  of  anti-Stokes  power  and  backward- radi¬ 
ated  Stokes  and  anti-Stokes  light.  Excessively  high  gain  and  the  abrupt 
threshold  of  SRS  were  key  difficulties  throughout  the  investigation. 

In  this  report  we  show  that  some  of  these  anomalies  (such  as 
the  anomalous  gain)  can  be  understood  to  be  a  result  of  light- trapping, 
as  we  had  previously  predicted.  1  The  presence  of  large  amounts  of 
anti-Stokes  power,  forward- to-backward  Stokes  ratio,  and  wrong  cone 
angles,  however,  remain  poorly  understood.  More  puzzling  in  some 
respects  is  the  behavior  of  SRS  in  some  gases  when  the  data  "suggest" 
that  beam-trapping  may  be  taking  place;  this  would  certainly  be  sur¬ 
prising  if  it  were  true.  From  our  theoretical  studies,  we  also  con¬ 
clude  that  any  liquid,  even  of  symmetric  molecules,  will  be  expected 
to  trap  when  the  total  power  in  the  material  (Raman  or  laser  pump) 
is  in  excess  of  10°  to  10?  W.  However,  the  use  of  gases  may  permit 
us  to  handle  much  more  power.  In  this  report  we  review  our  work  on 
SRS  and  the  beam- trapping  effects  which  have  now  become  almost  in¬ 
distinguishable  from  SRS,  although  there  are  clearly  two  different 
effects.  The  discussion  will  be  organized  with  the  experiment  and 
theory  separate,  although  the  two  could  not  be  and  were  not  separated 
at  the  time. 


1 


II. 


MODE  SELECTION 


At  the  beginning  of  this  contract,  considerable  emphasis  was 
placed  upon  the  fact  that  almost  all  of  the  theories  of  SRS  had  run  into 
severe  conflict  with  experiments.  One  of  the  main  reasons  for  this 
sudden  divergence  of  experiment  and  theory  was  the  possibility  of 
experiments  using  single  mode  GP  lasers  to  study  the  nonlinear 
effects,  as  described  in  Appendix  I.  All  of  the  early  nonlinear  work 
had  been  done  w’th  multimode  lasers  and  the  mode  structure  of  the 
laser  had  been  shown,  by  us  and  others,  to  cause  a  considerable 
spectral  broadening  of  the  anti-Stokes  radiation.  There  was  a  large 
anomaly,  not  generally  recognized,  in  the  Raman  gain;  the  first  ex¬ 
planation  used  the  multimode  arguments  that  had  been  used  before  in 
a  different  context.  The  multimode  theory,  of  course,  has  been  shown 
to  be  correct  as  far  as  it  goes,  but  our  single  mode  experience  showed 
that  the  theory  did  not  explain  the  magnitude  of  the  anomaly  in  the 
Raman  gain,  as  is  shown  in  Appendix  II  and  elsewhere.  2 

This  early  work  left  the  gain  anomaly  unresolved.  Fortu 
nately,  the  new  idea  of  light  trapping  proved  powerful  enough  to  pro¬ 
vide  a  resolution  ,  many  of  the  features  of  SRS  are  determined  by 
light- trapping  effects,  and  the  study  of  SRS  has  turned  almost  com¬ 
pletely  into  the  study  of  light- trapping. 
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m. 


STUDIES  OF  BEAM  TRAPPING 


Seif-focusing  and  trapping  phenomena  in  high  power  electro¬ 
magnetic  beams  were  foreseen  by  Askan'yan,  ^  Talanov,  4  and 
Chiao,  et  al .  ®  We  later  suggested  ^  that  this  mechanism  could  ex¬ 
plain  some  of  the  anomalies  observed  in  SRS,  especially  the  low 
thresholds,  and  the  effect  has  since  been  observed  in  many  liquids 
We  discuss  our  experimental  investigations  of  beam  trapping  in 
liquids,  describing  ( 1 )  the  time-integrated  pnuiographic  observations 
of  filaments  in  nitrobenzene  and  carbon  disulfide,  (2)  the  time-resolved 
photographic  studies  of  beam  trapping  in  nitrobenzene,  and  (3)  the 
polarization  effects  observed  in  beam  trapping  in  several  liquids.  This 
report  on  our  experimental  work  is  followed  by  a  discussion  of  how 
well  the  phenomenon  of  beam  trapping  explains  the  observed  anomalies 
in  SRS. 


A.  Time-Integrated  Photographic  Observations  of  Filaments 

The  experimental  apparatus  used  for  these  investigations  is 
shown  schematically  in  Figs.  1  and  2.  Figure  1  shows  the  mode-con¬ 
trolled  giant  pulse  ruby  laser  system  which  was  discussed  in  detail 
in  Appendix  I.  Figure  2,  with  the  STL  image  converter  camera  re¬ 
placed  by  an  ordinary  camera  and  the  0.  5  mm  slit  removed,  shows 
the  arrangement  we  used  for  observing  filaments.  In  this  arrangement 
the  beam  from  the  giant  pulse  laser  is  first  accurately  defined  by 
passage  through  a  1.  5  mm  diameter  aperture;  it  then  travels  through 
a  10  cm  long  sample  cell  where  SRS  occurs,  and  finally  goes  through 
the  following  optical  system,  along  with  the  stimulated  Raman  emission. 

The  optical  system  produces  dispersed,  magnified  images  of 
the  exit  end  of  the  sample  cell  on  the  film  plane  of  the  camera.  Photo¬ 
graphic  recording  of  these  images  allows  study  of  the  size,  number, 
contrast,  etc.  ,  of  filaments  at  the  laser  and  at  the  first,  second,  and 
(for  CSz)  third  order  Stokes  wavelengths.  The  wavelength  dispersion 
is  accomplished  by  imaging  the  end  of  the  sample  cell  through  a  small, 
high  speed,  diffraction  grating  spectrograph.  The  end  of  the  cell  is 
imaged  in  the  entrance  aperture  of  the  spectrograph,  the  dispersion  of 
the  latter  separates  the  images  for  the  various  wavelengths  ’n  the  exit 
plane.  The  camera  then  images  the  exit  plane  of  the  spectrograph 
onto  its  film  plane.  The  over-all  magnification  of  the  final  image  at 
the  film  plane  is  the  product  of  the  magnifications  of  the  first  two 
lenses  —  that  through  the  spectrograph  and  that  of  the  camera  —  and 
was  typically  eight  times  in  our  experiments.  The  resolution  obtained 
was  about  10  p  when  I-N  photographic  plates  were  used,  and  considera¬ 
bly  worse  for  infrared  Polaroid  film.  The  resolution,  cf  course, 
directly  determines  the  minimum  filament  size  which  can  be  inferred 
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Fig.  1.  Schematic  diagram  of  the  giant  pulse  laser  system. 
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Fig.  2.  Schematic  diagram  of  the  apparatus  used  for  observing  laser  and  Stokes 
filaments. 


from  these  observations.  Since  only  infrared  Polaroid  film  was  used 
in  the  experiments  with  CS2  we  do  not  have  definitive  data  for  that 
material,  and  will  proceed  to  discuss  the  more  complete  observa¬ 
tions  on  nitrobenzene.  There  is  no  indication  of  any  gross  difference 
for  CS2. 

The  observations  of  nitrobenzene  filaments  were  made  at  and 
somewhat  above  SRS  threshold.  Typically  about  10%  of  the  laser  pump 
beam  was  observed  to  be  trapped  into  small  filaments  about  30  to  50  fx 
ir  diameter.  Several  filaments  (5  to  10)  were  usually  observed  for  a 
single  firing  of  the  GPL  and  these  were  of  roughly  equal  intensity. 

The  integrated  intensity  of  these  laser  filaments  was  only  about  two  to 
three  times  that  of  the  inhomogeneous  background,  thus  making  obser¬ 
vation  of  them  somewhat  difficult.  Radiation  at  the  first  and  higher 
order  Stokes  wavelengths  occurred  only  in  filaments  which  were  usual¬ 
ly  smaller  than  the  laser  filaments.  Sizes  of  the  Stokes  filaments 
ranged  from  the  resolution  limit,  i.  e.  ,  about  10  fx,  up  to  about  30  |x 
in  diameter.  The  correlation  between  the  laser  and  Stokes  filaments 
was  very  good,  indicating  that  SRS  was  indeed  occurring  in  the  trapped 
laser  filaments  and  that  the  Stokes  radiation  thus  produced  was  also 
trapped  in  the  same  filamentary  region.  Occasionally,  a  Stokes  fila¬ 
ment  was  observed  without  a  corresponding  laser  filament.  This  can 
be  attributed  to  saturation  of  the  laser  pump  radiation  in  that  filament, 
with  almost  total  conversion  to  Stokes  radiation. 

The  number  of  laser  and  Stokes  filaments  is  observed  to  in¬ 
crease  uniformly  with  increasing  pump  power  above  threshold.  The 
power  in  a  single  filament,  however,  appears  to  be  nearly  a  constant. 


B. 


Time-Resolved  Photographic  Studies  of  Filaments 


The  phenomenon  of  beam  trapping  can  explain  the  anomalously 
low  thresholds  for  SRS  if  the  intensity  in  the  trapped  filament  is 
sufficiently  higher  than  the  background  intensity.  As  noted  above,  we 
observed  a  time-integrated  enhancement  of  intensity  of  only  two  to 
three  times  in  the  filaments,  whereas  approximately  a  factor  of  100 
is  needed  to  explain  the  discrepancy  with  theory.  *  The  absence  of 
such  large  enhancements  in  the  above  data  can  be  explained  if  a  typical 
filament  is  present  for  a  time  short  compared  to  the  total  pump  pulse 
length.  Since  the  latter  is  typically  5C  nsec,  the  filament  duration 
should  be  only  one  or  a  few  nanoseconds  to  explain  the  observed  thres¬ 
holds.  In  order  to  check  this  hypothesis,  we  undertook  an  investi¬ 
gation  of  the  time-resolved  laser  and  Stokes  emission. 


The  apparatus  for  these  experiments, shown  schematically  in 
Fig.  2,  is  the  same  as  that  for  the  time -integrated  observations  des- 
cr'bed  above  except  for  the  STL  image  converter  camera  and  a  0.  5  mm 


slit  at  the  entrance  aperture  of  the  spectrograph.  In  order  to  get 
maximum  time  resolution,  the  fastest  streak  mode  of  operation  of 
the  STL  camera  was  used.  This  was  a  50  nsec  sweep  and  resulted 
in  a  time  resolution  of  about  1  nsec.  In  some  experiments  the 
0.  5  mm  slit  was  omitted,  s~>  that  a  smear  photograph  was  obtained 
instead  of  a  streak  photograph. 

The  short  streak  duration  and  the  resulting  high  writing  rate 
led  to  several  experimental  difficulties  which  limited  the  amount  and 
qualify  of  the  data.  The  experiments  required  operation  of  the  STL 
camera  at  the  limits  of  its  capabilities.  For  example,  some  com¬ 
promise  had  to  be  made  between  adequate  exposure  and  overloading 
the  electron  optics  in  the  STL  camera.  Two  difficulties  arising  from 
this  were  rather  poor  resolution  compared  with  that  of  the  static 
photographs,  and  3ome  ambiguity  in  interpretation  of  the  observed 
data,  i.  e.  ,  in  the  decision  as  to  whether  some  feature  was  real  or 
caused  by  distortion  of  the  electron  optics.  By  taking  some  data  at 
very  low  exposures  and  using  the  static  photographs  as  a  guide,  we 
feel  we  have  generally  been  able  to  make  good  interpretations  of 
the  observed  data.  Another  major  difficulty  arising  from  the  short 
streak  duration  was  in  proper  triggering  of  the  STL  camera.  The 
usual  trigger  circuits  were  found  to  be  too  slow,  and  various  other 
triggering  schemes  were  used.  The  most  useful  is  that  shown  in 
Fig.  where  the  output  of  a  planar  ("bomb")  photodiode,  driven 
very  hard  by  part  of  the  GPL  output,  is  vpplied  directly  to  trigger 
the  camera.  With  this  arrangement  we  were  usually  able  to  observe 
the  significant  parts  of  the  laser  pump  pulse. 

Our  time-resolved  observations  were  on  the  laser  and  first 
Stokes  radiation  from  nitrobenzene.  As  in  the  above  observations, 
the  laser  beam  was  linearly  polarized.  The  laser  filaments  appeared 
at  different  times  lasting  typically  from  3  to  6  nsec.  The  first  Stokes 
filaments  were  observed  to  persist  for  a  shorter  time  of  about  1  to  3 
nsec.  The  peak  intensity  of  the  laser  filaments  relative  to  the  back¬ 
ground  intensity  could  not  be  accurately  determined  because  of  the 
small  contrast  range  inherent  in  the  experiments,  but  was  at  least 
one  order  of  magnitude  higher.  Comparison  of  the  laser  and  Stokes 
filament  durations  and  the  laser  pulse  duration  indicates  that  these 
observations  explain  the  threshold  discrepancy  at  least  to  within  a 
factor  of  two. 

In  the  time-resolved  photographs  we  have  again  occasionally 
observed  Stokes  filaments  without  corresponding  laser  filaments, 
and  very  bright  laser  filaments  with  no  evident  associated  Stokes 
emission.  In  addition,  there  appears  to  be  some  spatial  "wandering" 
of  the  observed  filaments  that  is  not  a  result  of  distortion  in  the 
electron  optics.  Although  these  issues  are  not  definitely  settled,  we 
feel  the  effects  are  probably  due  to  saturation  of  the  laser  pump  and 
the  short  time  scale  dynamics  of  filament  formation. 


c. 


Polarization  Effects  in  Beam  Trapping 


The  physical  mechanism  which  has  usually  been  accepted  as 
responsible  for  beam  trapping  in  liquids  is  molecular  reorientation 
in  the  strong  optical  fields;  this  is  generally  teamed  the  ac  Kerr 
effect.  Comparisons  of  the  Kerr  constants  with  observed  thresholds 
in  liquids  which  readily  exhibit  SRS  have  been  usually  favorable. 

We  have  attempted  to  verify  this  molecular  reorientation  hypothesis 
by  observing  the  filaments  induced  by  pump  light  which  is  linearly 
or  circularly  polarized,  and  by  seeing  whether  the  observed  thres¬ 
holds  for  these  two  cases  behave  according  to  the  theory  of  reorien¬ 
tation.  We  have  made  photographic  and  photometric  studies  of  the 
polarization  dependence  of  beam  trapping  in  several  liquids.  We 
have  found  an  instability  in  the  nonlinear  propagation  of  circularly 
polarized  light  which  has  made  interpretation  of  the  threshold 
measurements  difficult.  However,  the  origin  of  this  instability  and 
its  consequences  are  interesting.  A  paper,  "The  Self-Focusing  of 
Light  of  Different  Polarizations,  "  which  discusses  these  results  and 
which  will  be  published  in  the  Fourth  Quantum  Electronics  Conference 
Proceedings,  is  reproduced  here  as  Appendix  III. 

1.  Photographic  Techniques 

The  apparatus  for  the  photographic  experiments  is 
shown  schematically  in  Fig.  3,  and  is  the  same  as  that  used  for  the 
static  photographic  studies,  except  that  a  quarter-wave  plate  to  fix 
the  polarization  of  the  laser  pump  beam  and  a  Fresnel  rhomb  and 
Wollaston  prism  to  analyze  the  polarization  of  the  radiation  exiting 
from  the  sample  cell  have  been  added.  The  camera  thus  records 
dispersed  images  of  the  cell  end,  with  polarization  information  at 
each  wavelength.  By  setting  the  Wollaston  prism  properly,  the 
analyzer  system  gives  either  the  components  of  the  radiation  in  two 
orthogonal  linear  polarizations  or  the  components  in  the  right  and 
left  circular  polarizations. 


Photographic  observations  were  made  on  nitrobenzene  and  CS^. 
If  the  imput  beam  was  linearly  polarized,  the  laser  background,  laser 
filaments,  and  Stokes  filaments  were  all  like -polarized  to  within  one 
part  in  10  ,  the  limit  of  our  instrumentation.  However,  if  the  input 
beam  was  circularly  polarized,  each  Stokes  order  for  nitrobenzene 
contained  equal  amounts  of  both  circular  polarizations  over  the  entire 
range  of  observation  of  about  three  orders  of  magnitude  in  the  total 
first  Stokes  intensity.  Furthermore,  the  right  and  left  circularly 
polarized  images  at  any  one  Stokes  frequency  were  identical  in  every 
detail  of  intensity  and  pattern.  The  laser  background  appeared  en¬ 
tirely  in  the  input  polarization,  while  the  laser  filaments  appeared 
in  both  circular  polarizations.  The  laser  filaments  also  seemed  to 
be  equally  intense  and  otherwise  identical  in  both  circular  polari¬ 
zations,  although  this  was  difficult  to  ascertain  with  the  laser  back- 
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beam  trapping. 


ground  appearing  in  the  one  circular  polarization.  When  the  Wollaston 
prism  was  rotated  to  test  for  orthogonal  linearly  polarized  components 
with  the  same  (circular)  input  polarization,  the  nitrobenzene  Stokes 
filament  patterns  for  the  two  linear  polarizations  at  one  frequency  did 
differ  in  detail,  with  variations  in  intensity  between  corresponding 
filaments  and  many  cases  of  unpaired  filaments  in  one  or  the  other  of 
the  polarizations.  As  expected,  the  laser  background  appeared 
equally  strong  in  both  linear  polarizations,  making  analysis  of  the 
laser  filaments  difficult.  The  total  Stokes  intensity  was  about  equal 
in  each  linear  polarization  at  one  frequency.  These  observations  led 
to  the  interpretation  that  the  circular  input  beam  produces  filaments 
in  nitrobenzene  which  are  linearly  polarized,  with  a  random  direction 
of  polarization  which  varies  from  filament  to  filament.  The  theoretical 
basis  for  these  observations  is  discussed  below. 

Similar  experiments  with  circularly  polarized  pump  light  in 
CS2  showed  some  slightly  different  results.  There  was  incomplete 
depolarization,  with  about  a  third  of  the  filament  energy  appearing 
in  the  opposite  circular  polarization  to  that  of  the  laser  and  first 
Stokes  filaments.  The  partition  of  energy  was  more  nearly  equal  :n 
the  second  Stokes  output  and  became  essentially  equal  in  the  third 
Stokes  output.  This  can  be  taken  as  evidence  that  the  various  Stokes 
orders  are  iteratively  generated,  i.  e. ,  that  the  laser  produces  the 
first  Stokes,  which  in  turn  produces  the  second  Stokes,  etc. 

Since  a  wide  range  of  Stokes  intensity  could  not  be  covered  in 
the  photographic  observations,  estimates  of  the  relative  thresholds 
for  circular  and  linear  input  polarizations  are  not  very  accurate.  By 
counting  the  number  of  observed  Stokes  filaments  as  a  function  of 
laser  power,  we  found  an  approximate  ratio  of  circular  threshold  to 
linear  threshold  of  about  3.  8  for  nitrobenzene  and  2.  5  for  CS>2»  which 
can  be  compared  with  the  theoretical  value  of  4  for  the  ac  Keir  effect  and 
with  the  more  accurate  photometric  results  below. 

2.  Photometric  Techniques 

For  the  photometric  measurements,  the  unfocused 
beam  of  a  giant  pulse  ruby  laser  was  first  passed  through  crossed 
Gian  polarizers  which  served  as  a  variable  attenuator  and  supplied  a 
purely  linearly  polarized  output  beam.  The  unfocused  beam  then 
passed  through  a  quarter- wave  plate  and  a  10  cm  cell  containing  the 
liquid  in  which  self-trapping  was  to  be  studied.  A  portion  of  the  output 
beam  from  this  cell  was  split  off  and  monitored  by  silicon  photodiodes 
for  Stokes  energy  content  after  passing  through  a  Corning  5-56  filter, 
two  Wratten  89B  filters,  and  appropriate  attenuators.  When  the  liquid 
cell  was  removed,  the  laser  produced  no  signal  at  the  photodetector, 
thus  assuring  that  only  Stokes  light  was  being  measured. 


12 


Typical  results  for  the  Stokes  output  powers  plotted  as  a 
function  of  either  linearly  or  circularly  polarized  inputs  are  given  for 
benzene,  nitrobenzene,  and  CS^  in  Fig.  4.  From  such  data  the  trap- 
ping  threshold  ratios{circular  to  linear)  for  various  liquids  are  given 
in  Table  I  of  Appendix  III.  These  ratios  are  seen  to  be  consistently 
well  below  the  value  4  expected  if  circularly  polarized  light  would 
produce  circularly  polarized  filaments  via  the  ac  Kerr  effect. 

The  depolarization  of  the  Stokes  Raman  output  from  nitroben¬ 
zene  was  also  studied  photometrically  by  using  the  polarization  analysis 
procedure  of  Fig.  3.  The  exit  window  of  the  cell  was  imaged  through 
a  Fresnel  rhomb  and  a  Wollaston  prism  to  form  two  images  on  a 
boron  nitride  surface.  The  images  were  monitored  separately  by 
silicon  photodiodes  through  filters.  Depending  on  the  relative  orien¬ 
tation  of  the  Wollaston  prism  to  the  Fresnel  rhomb,  signals  produced 
by  the  two  images  were  proportional  either  to  the  right  and  left  circu¬ 
lar  component  intensities  or  to  the  intensities  of  the  components 
linearly  polarized  parallel  and  perpendicular  to  the  original  laser 
polarization.  When  the  cell  was  removed,  the  signals  showed  that  only 
the  nominal  input  polarization  (either  linear  or  circular)  was  present 
to  within  one  part  in  200,  the  limit  of  resolution.  When  Stokes  light 
was  produced  by  linearly  polarized  input,  it  was  found  to  be  linearly 
polarized  to  within  one  part  in  200;  when  it  was  produced  by  a  circularly 
polarized  input,  it  was  found  to  contain  equal  intensities  of  both  circular 
polarizations  and  of  both  linear  polarizations.  These  results  applied 
over  the  entire  power  range  where  Stokes  light  could  be  observed 
(roughly  seven  orders  of  magnitude  variation  in  Stokes  output),  and 
are  consistent  with  the  photographic  observations.  Also,  apertures 
ranging  from  l/2  to  5  mm  in  diameter  placed  in  the  laser  beam  did 
not  affect  these  findings,  even  when  the  cell  containing  the  liquid  was 
placed  in  the  far  field  of  the  1/2  mm  aperture. 


D.  Resolution  of  Anomalies  in  SRS 


The  observation  of  the  self-trapping  now  associated  with  most 
SRS  in  liquids  was  a  major  step  forward  in  the  understanding  of  the 
SRS  process  in  these  media.  We  now  discuss  how  well  the  concept  of 
self-trapping  can  explain  the  anomalies  which  have  been  observed.  The 
general  conclusion  is  that  the  phenomenon  satisfactorily  explains  some 
features,  explains  some  other  features  in  principle  but  not  in  detail, 
and  leaves  some  questions  still  unanswered,  even  conceptually. 

The  most  striking  of  the  puzzling  anomalies  are  the  low  thres¬ 
holds  and  the  large  wavelength  spreads  with  multimode  pumping.  These 
are  also  the  features  which  have  received  the  best  explanations.  The 
time-resolved  observations  of  filaments  in  nitrobenzene  discussed 
above  showed  that  the  low  thresholds  can  be  quantitatively  explained 
as  resulting  from  transient  beam  trapping,  Bloembergen  and  Lallemand^ 
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INCIDENT  ENERGY 


Fig.  4.  Plots  of  the  Raman  Stokes  output  energy  versus 
ruby  laser  energy  incident  on  a  10  cm  cell  of  (a) 
nitrobenzene,  (b)  carbon  disulphide,  and  (c) 
benzene.  Triangular  points  represent  linearly 
polarized  incident  light,  circular  points  represent 
circularly  polarized  light.  The  zero  of  the  log 
Stokes  output  energy  corresponds  roughly  to  a 
conversion  of  10“^  of  the  laser  energy  to  Stokes 
energy.  The  incident  energy  was  not  calibrated 
absolutely.  The  dotted  lines  represent  pulses 
producing  no  detectable  Stokes  output. 


14 


have  shown  that  frequency  broadening  can  be  explained  in  terms  of 
orientation  effects  in  anisotropic  molecules.  These  results  were 
sufficiently  qualitative  that  another,  similar  process  may  be  found 
pertinent,  just  as  we  feel  that  orientation  effects  may  not  completely 
explain  the  beam  trapping  phenomenon.  In  any  case,  the  frequency 
spreading  is  definitely  related  to  the  beam  trapping. 

The  forward  Stokes  to  backward  Stokes  ratio  measurement  is 
generally  complicated  by  feedback  to  the  ruby  laser  via  stimulated 
Brillouin  scattering  (SBS).  Also,  the  coupling  to  the  anti -Stokes 
radiation  is  very  important  to  the  development  of  the  Stokes  radiation, 
and  the  former  is  not  well  understood.  ^  Therefore,  we  must  conclude 
that  the  forward-backward  Stokes  ratio  is  not  completely  understood, 
but  we  may  hope  it  can  be  understood  within  present  concepts.  Essen¬ 
tially  the  same  comment  can  be  made  regarding  the  importance  of 
regeneration,  the  Stokes  to  anti-Stokes  ratio,  and  the  anti-Stokes 
emission  angles.  While  the  single-pass  gain  can  be  sufficiently  large 
under  trapped  conditions  to  obtain  the  observed  output  from  a  single 
pass,  the  details  of  the  frequency  composition  and  coherence  pro¬ 
perties  of  the  SRS-induced  radiation  are  not  well  understood.  Like¬ 
wise,  the  anomalously  large  amount  of  anti-Stokes  radiation  produced 
and  the  "wrong"  anti-Stokes  emission  angles  have  been  shown  by 
Garmire1^  and  by  Bret  and  Denariez1*  to  be  intimately  associated 
with  beam  trapping,  but  the  processes  invovled  are  not  at  all  well 
understood. 

Probably  the  most  pressing  remaining  questions  concern 
SRS  in  gases.  The  observation  by  Terhune1*  of  large  amounts  of 
anti-Stokes  radiation  in  hydrogen  gas,  and  the  measurements  by 
Mayer15  of  SRS  conversion  properties  in  hydrogen  very  similar  to 
those  observed  in  liquids  showing  beam  trapping  are  extremely  puzzl¬ 
ing.  There  is  no  known  process  which  could  furnish  a  nonlinear  index 
large  enough  to  obtain  beam  trapping  in  these  cases,  and,  as  far  as 
is  known,  nothing  but  beam  trapping  has  the  potential  for  explaining 
these  effects.  Thus,  many  questions  remain  open,  and  it  appears 
that  the  investigation  of  gases  and  their  comparison  with  liquids  will 
be  a  very  fruitful  avenue  of  approach  to  understanding  more  about  the 
complex  SRS  process. 

E.  Theory  of  Beam  Trapping 

1.  Introduction 


The  "self-focusing"  or  "self-trapping"  of  the  pump  light 
or  the  Raman  scattered  light  in  a  Raman  laser  tends  to  destroy  the 
coherence  of  the  Raman  output.  Therefore,  self-focusing  effects  mustbe 
prevented  if  a  high  coherence  high  energy  Raman  output  is  to  be  obtained. 
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It  is  generally  felt  that  index  changes  from  electros triction  and  heating 
which  involve  macroscopic  mass  transport  in  the  Raman  media  will 
not  develop  significantly  in  times  short  compared  with  10"?  -  10"8  sec. 
However,  much  faster-acting  effects  which  do  not  involve  macroscopic 
mass  transport  may  cause  index  changes  and  self-trapping.  The  most 
studied  of  the  fast  self-trapping  effects  is  that  which  arises  from  the 
reorientation  of  asymmetric  molecules  in  the  strong  optical  fields  and 
is  called  the  ac  Kerr  effect.  The  thresholds  for  self- trapping  that  are 
expected  theoretically  from  the  ac  Kerr  effect  are  in  rough  quantita¬ 
tive  agreement  with  the  observed  thresholds  in  such  asymmetric  mole¬ 
cular  liquids  as  CS2*  nitrobenzene,  benzene,  toluene,  and  benzoyl- 
chloride.  These  power  thresholds  lie  in  the  10  to  100  kW  region,  or, 
in  terms  of  pulse  energies,  in  the  10“^  -  10-4  J  region.  These  thres¬ 
hold  energy  values  give  essentially  an  upper  limit  on  the  Raman  output 
power  that  can  be  obtained  in  an  undistorted  wave  train.  Therefore, 
for  the  purposes  of  the  present  study,  we  would  like  to  see  if  some 
Raman  active  liquids  could  be  found  which  have  three  or  more  orders 
of  magnitude  higher  trapping  thresholds.  One  might  hope  that  a 
perfectly  symmetric  molecule  which  can  exhibit  no  ac  Kerr  effect 
might  therefore  have  a  trapping  threshold  (in  liquid)  many  orders  of 
magnitude  greater  than  those  liquids  mentioned  above.  In  this  section 
we  are  concerned  with  the  theoretical  investigation  of  this  possibility. 
The  main  result  is  that  no  such  hope  exists;  that  is,  no  Raman  active 
liquid  near  room  temperature  can  be  expected  to  sustain  a  wave  of 
around  10°  W  or  greater  without  deterioration  of  beam  coherence  due 
to  self-focusing  effects. 

If  the  main  contribution  to  the  fast- responding  part  of  the  non¬ 
linear  index  of  a  liquid  of  symmetric  molecules  were  a  result  of  the 
nonlinear  polarizability,  perhaps  an  increase  in  trapping  threshold 
powers  of  three  or  more  orders  of  magnitude  (10?  to  10*0  W)  could  be 
expected.  However,  there  is  another  contribution  to  the  nonlinear 
index  which  is  much  larger,  being  of  the  order  of  the  reorientation 
contribution  for  asymmetric  molecules,  and  remaining  large  for  sym¬ 
metric  molecules.  This  contribution  is  from  the  local  molecular  re¬ 
distribution  in  space  that  must  occur  in  the  presence  of  strong  optical 
fields. 

The  probability  of  finding  a  certain  spatial  configuration  of 
molecules  in  a  liquid  depends  on  the  energy  of  that  configuration.  When 
a  large  optical  field  is  present  in  the  liquid,  the  relative  energies  of 
various  configurations  will  change  owing  to  the  added  dipole-dipole 
interactions.  Configurations  having  the  largest  polarization  become 
energetically  more  favorable  and  hence  more  probable,  and  an  increase 
in  the  index  of  refraction,  similar  to  that  from  reorientation,  must 
result  from  the  spatial  redistribution.  This  can  happen  either  with  a 
macroscopic  density  change  (electrostriction)  and  mass  transport  or 
without  either.  Here  we  compute  the  nonlinear  index  change  from 
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microscopic  redistributions  which  involve  no  macroscopic  density 
changes  in  the  liquids  and  hence  can  occur  in  times  very  short 
compared  with  a  nanosecond  pump  pulse  duration.  It  is  such  spatial 
redistribution  which  will  determine  the  short-pulse  self- trapping 
thresholds  of  fluids  of  symmetric  molecules  rather  than  electronic 
or  eiectrostrictive  effects.  In  fact,  we  will  find  that  even  for  such 
asymmetric  molecules  as  nitrobenzene  and  benzene  redistribution 
effects  are  as  important  as  the  ac  Kerr  effect  in  self-trapping. 
Examination  of  the  predicted  relative  thresholds  of  all  liquids  for 
which  measuremei  have  been  made  at  Hughes  or  reported  else¬ 
where  show  all  behaving  in  accord  with  the  estimates  made  below 
from  combined  reorientation  and  redistribution.  For  example,  the 
symmetric  molecule  CCI4  is  predicted  to  have  a  lower  threshold 
than  the  asymmetric  molecules  chloroform,  acetone,  acetic  acid, 
and  water,  and  a  higher  threshold  than  benzene,  nitrobenzene,  and 
toluene  as  is  observed.  It  is  further  suggested  that  some  other 
symmetric  molecules  such  as  tetraisopropyl  lead  which  have  not 
yet  been  studied  experimentally  will  also  self-trap  light  at  thresholds 
as  low  as  those  at  some  asymmetric  molecules  which  have  exhibited 
the  effect. 

In  the  course  of  estimating  the  nonlinear  index  of  liquids 
of  symmetric  molecules  we  derive  a  rigorous  lower  bound  for  the 
linear  index  within  the  framework  of  classical  statistical  mechanics. 
This  quite  strong  inequality  is  of  interest  because  it  contradicts 
(by  a  very  small  amount)  recent  experimental  results  on  argon, 
methane,  CF^,  and  other  symmetric  materials. 

2.  Formulation 


We  calculate  here  the  nonlinear  dielectric  constant 
as  if  the  fields  were  static.  The  optical  dielectric  constant  can  be 
obtained  from  the  static  result  simply  by  replacing  the  static  polar¬ 
izability  by  the  optical  polarizability  since  permanent  moments  have 
no  effect  on  the  optical  properties.  Also,  die  square  of  the  static 
electric  field  must  be  replaced  by  the  average  over  many  cycles  of 
the  square  of  the  optical  field.  That  the  intermolecular  forces  and 
molecular  torques  are  essentially  the  same  for  optical  as  for  static 
fields  is  because  the  motions  of  the  molecules  in  the  liquid  cannot 
respond  to  the  forces  existing  at  optical  frequencies,  but  only  to  the 
relatively  low  frequency  components  of  force.  The  logic  of  die  ex¬ 
tension  of  the  static  results  to  the  optical  case  will  become  clear  in 
the  following  development.  At  the  cost  of  considerable  complexity, 
the  development  can  in  fact  be  made  mathematically  rigorous  for 
optical  fields  throughout  and  achieve  essentially  the  same  result. 
However,  the  simpler  treatment  that  is  available  for  the  nonlinear 
static  dielectric  constant  will  be  pursued  here  because  of  its  greater 
physical  clarity. 
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To  find  the  static  dielectric  constant,  we  compute  the  free 
energy  3*  of  an  ellipsoidal  sample  placed  in  a  uniform  x-directed 
electric  field  E0  such  that  an  axis  of  the  ellipse  is  parallel  to  the 
field.  Then  the  total  electric  dipole  moment  of  the  sample  M  is 
obtained  from  the  well  known  relation  (proved  later) 

M  =  -  ds3»/dE0  .  (1) 

Because  the  sample  is  ellipsoidal  and  the  external  field  uniform,  the 
polarization  per  unit  volume  P  inside  the  sample  will  be  uni  orm 
and  x-directed.  We  assume  the  sample  to  have  a  fixed  volume  V. 
Therefore, 

P  =  M/V  .  (2) 

We  desire  the  nonlinear  index  only  to  second  order  in  the  electric 
fields  so  we  compute  3*  to  fourth  order  in  the  fields  and  obtain  a 
result  of  the  form 

31  =  (bE|/2  +  c  E|/4)V  (3) 

from  which  we  have 

P  =  b  E0  +  c  E*  .  (4) 

In  an  ellipsoidal  sample  uniformly  polarized  along  a  principal 
axis,  the  relation  of  the  polarization  density  to  the  macroscopic 
electric  field  E  is  generally  written^ 

E  =  E  -  LP  (5) 

o  '  ' 

where  L  is  called  the  depolarization  factor  for  that  axis  of  the 
ellipsoid.  The  factor  L  can  have  values  between  0  and  4tt  depend¬ 
ing  only  on  the  shape  of  the  ellipsoid.  ^  If  we  define  a  linear  sus¬ 
ceptibility  X  and  a  nonlinear  susceptibility  t|  by  the  relation 

P  =  XE  +  t|E^,  (6) 

a  comparison  of  (4)  and  (6)  with  the  aid  of  (5)  shows  that 

X  =  b/(l  -  Lb)  (7) 

and 

r|  -  c(l  -  Lb)-4  =  c(  1  +  LX)4  .  (8) 
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The  nonlinear  dielectric  constant  cXTT  is 

N  L 


«NL  =  1  +  4tt  P/E  =  1  +  4tt  (X  +  t|  E4), 


(9) 


and,  in  extrapolating  our  results  to  the  optical  case,  we  would  there- 
fore  have  an  index  of  refraction  n^^  given,  to  second  order  in  E,  by 


nNL  =  n  +  (2Tr/n)  tVE 


(10) 


where  E  is  the  rms  electric  field  averaged  over  many  optical  cycles 
and  n  is  the  linear  index  of  refraction.  Following  the  standard 
practice,  we  define  a  "nonlinear  index"  n^  by 


n. 


5  2irr)/n  . 


(ID 


The  threshold  power  required  to  begin  trapping  a  smooth  beam 
is5  c 


P  a  (1.  22  X.  )4  c  /( 1 28  nj 
c  o  o'  2 


(12) 


where  \Q  and  cc  are  the  vacuum  wavelength  and  veloc  iy  of  light 
respectively.  We  now  have  all  the  necessary  relations  to  connect 
the  nonlinear  index  effects  with  the  electric  free  energy^. 

In  computing  the  free  energy we  assume  that  we  are  dealing 
with  a  classical  fluid  and  use  classical  statistical  mechanics.  We  must 
therefore  formulate  the  total  potential  energy  v  of  the  N  molecules  in 
the  sample  when  they  are  placed  at  positions  ,  ...  &N,  have  electric 

dipole  moments  ml,  . . .  m^,  and  are  oriented  at  the  Euler  angles 
which  we  symbolize  by  Qj,  ...  Qpq. 

We  anticipate  that  the  electronic  nonlinearities  will  be  small 
compared  to  those  with  which  we  are  concerned,  and  assume  that  the 
internal  energy  v-y  of  the  y^1  molecule  is  quadratic  in  the  components 
of  its  electric  dipole  moment  m\  If  mY  are  the  components  of  this 
moment  along  the  principal  axes  a  =  l,^,  3,  of  its  polarizability 
ellipsoid,  then  17 


\  ■  7  i  =  i  <mI>2  /“a  (13> 

where  aj,  a^,  <13  are  the  polarizabilities  of  the  molecule  along  its 
principal  axes.  The  components  m^  are  related  to  the  components 
mY  measured  in  a  laboratory-fixed  coordinate  system  by* 


*We  will  consistently  use  the  subscripts  a  through  h  to  denote  the 
space  components  of  a  vector  or  matrix  referred  to  the  principal 
axes  1,  2,3  of  the  polarizability  ellipsoid  fixad  in  the  molecule. 
Subscripts  from  i  onward  always  denote  space  components  referred 
to  axes  fixed  in  the  laboratory. 
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(14) 


=  R.  (Q  )mY  . 
i  la  v  a 

Th>  reneated  space  indices  here  and  elsewhere  are  assumed  to  be 
summt  '  The  matrix  R^a  (Q^)  =  (^y  V  represents  the  rotation 

ope  rst  £c  or  the  Euler  angles  Q.y  Y^*1  molecule,  and  will  be 

abbreviated  RjY  . 

The  interaction  energy  v  ^  between  two  molecules  at  positions 
r^  and  r^  v  ill  be  assumed  to  De^ 


YU 


=  v  ...  +  m!  m^D^ 


oyn 


,V 

i 


J  ij 


(15a) 


where  vQ  „  is  a  hard  core  repulsive  potential  and  the  remaining 
term  is  tne  dipole-dipole  interaction  energy: 


(6. .  5,, 
dyh  _  m 
ij 


36ik  »«><*£ 


;>  <rVi  - 


V  M- 
r  *  -  r 


(15b) 


The  m*  may  be  taken  to  be  the  independent  harmonic  "internal" 
coordinates  olthe  molecules  which  describe  the  electronic  polarization. 
That  is,  it  is  appropriate  to  integrate  over  all  values  of  these  coordi¬ 
nates  in  evaluating  the  free  energy  from  the  potential  energy  function. 
The  free  energy  «^of  the  liquid  in  the  external  field  may  therefore  be 
written 


where 


sJJJ  dUf  di«  f  <Hmaf e_|3v  <16a) 


=  2(-hYE  +v)+i 
y  x  o  y  2 


2  v 
Y.H-  YM- 


(16b) 


is  the  total  potential  energy  of  a  configuration  of  the  positions,  moments, 
orientations  of  the  molecules.  J  djr^  symbolizes  the  integration 
Jdr^1  .  .  .  Jdr^*  over  the  entire  volume  V  of  the  ellipsoid  of  every 
molecular  position  and  similarly  for  all  molecular  orientations]  12  [  and 
internal  "nolecular  coordinates  |  ma  [  .  (3  is  defined  as  the  inverse  of 

(the  temperature  times  Boltzmann's  constant). 

This  completes  the  formulation  of  the  static  dipolar  free  energy 
suitable  for  classical  liquids.  The  relations  (1)  through  (12)  connect 
the  dependence  ofvT’on  an  externally  applied  electric  field  EQ  with  the 
threshold  power  for  self-focusing  or  self-trapping  which  we  wish  to 
predict.  It  remains  to  carry  through  the  evaluation  of  the  free  energy^ 
starting  from  the  basic  definition  (16). 
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3. 


A  Variational  Principle  for  Calculation 


_  In  order  to  obtain  from  (16)  an^ccurate  estimation 

of  \T  for  liquids  it  does  not  suffice  to  expand  in  powers  of  the 
number  density  p  =  N/V.  We  are  therefore  going  to  employ  a 
variational  approach  which  obtains  an  upper  bound  for  .  Param¬ 
eters  will  exist  in  the  upper  bound  which  can  be  varied  to  obtain 
a  lowest  upper  bound  within  any  approximation  scheme. 


To  obtain  a  variational  principle,  we  note  that  the  total  po¬ 
tential  energy  v  is  a  quadratic  function  of  the  dipole  moment  coor¬ 
dinates  m^.  This  means  that  if  we  define  new  coordinate  variables 
of  integration  by 


(17) 


v 

where  the  p'  are  those  values  which  minimize  v  for  a  given  positional 
and  orientational  configuration,  then  we  can  rewrite  (16)  as 


=  Voyp  ■(3v'  ‘  P 


min; 


>f] 


where  (see  footnote  on  p.  ) 


(18) 


t  E0> 


Uv 


+ A  2  \>y' 


p‘  E  +  £  P' 

x  o  2  a  a 

'W.K«KC. 


(19) 


is 


the  function  which  is  minimized  by  our  choice  of  the  p 


V 

a 


and 


exp[-(3v'(]r  )]  (v  a[  expt- PS(v^  /a^  + 

z  v»)  ] . 

p  1  u  J 


(20) 


The  symbol]  [  represents  the  set  of  all  values  of  the  pa.  The 
intermolecular  potential  v'  defined  by  (20)  cmnot  depend  on  the  electric 
field  E0  (or  the  ]  ji  }  would  not  minimize  ft)  and  represents,  of  course, 
the  long  range  van  der  Vaal's  potential  due  to  induced  dipole-induced 
dipole  interactions.  This  potential  has  been  analyzed  in  the  literature^ 
and  we  need  not  discuss  it  further  here. 
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It  is  too  difficult  to  find  an  exact  explicit  expression  for 
[  min  1^1  ] ,  but  if  we  insert  wrong  values  for  the  j  ja  J»  into  $  we 

at  leaiat  obtain  an  estimate  for  the  field-dependent  part  ,j£  of  >3  that 
is  too  high.  We  therefore  have  the  useful  inequality  that  for  any 
|  jx[  whatever 

v3e  <  F  (21) 

where 

e_f3F  5  JJdW  dinl  (22a> 

and 


e*p(-pswv0v|1-pv') _ 

|d|r}djn}  exp 


(22b) 


With  this  inequality  we  may  proceed  with  a  variational  calculation  of 
the  desired  field  dependent  partir^,  of  the  free  energy. 


The  Effective  Field  Approximation 


We  will  use  a  local  field  approximation  in  which  we  use 
the  approximate  values 


X. 

^a  = 


a  R  Y  E' 
a  xa 


(23a) 


or  equivalently 

\iy.  =  R.Y  q  R  Y  E'  (23b) 

j  ja  a  xa 

as  if  each  molecular  moment  saw  only  an  effective  x-directed  local 
field  of  magnitude  E\  We  will  then  adjust  E'  to  minimize  the  re¬ 
sulting  upper  bound  F  for  the  free  energy  (calculated  to  fourth  order 
in  E').  The  F  thus  obtained  is  a  function  of  the  shape  of  the  ellipsoidal 
sample,  but  the  dielectric  constant  derived  from  it  via  (9)  is  not.  (Other¬ 
wise  we  could  vary  the  shape  parameters  to  improve  further  the  accuracy 
of  the  estimate  for  the  dielectric  constant.)  Using  (23)  in  (19)  gives 


22 


where 


and 


e‘PF  =«e-fJ(u2  '  ul)» 

(24) 

“i  2  zv  S  •  a,Y  •  S'  • 

(25a) 

G  5  E  -  E'/2 

(25  b) 

is  E'  .  a Y  V  .  a V  .  E'  . 

2  y9V  ~  **  ***  **  rsj 

(26) 

We  have  introduced  the  usual  shorthand  notation  for  multiplication 
among  the  3x3  space  matrices 


la 


a  R.T 
a  ja 


A  A 

and  the  three- vectors  Eo  =  x  E0  and  jg'  =  x  E'.  The  double  average 
brackets  «  »  in  (24)  symbolize  the  integrations  over  and  j  ft  [ 

weighted  by  w( |x[)  which  is  normalized  so  that  «  1  =  1.  We 

will  often  use  only  the  single  average  bracket  when  either  the  angle  or 
position  average  is  required  alone;  which  average  is  intended  will  be 
clear  from  the  context.  In  either  case,  <  1  >  =  1. 

We  expect  (as  we  will  find)  that  the  best  local  field  E'  will  be 
of  the  order  of  E0.  Therefore,  eq.  (24)  implies  that,  to  fourth  order 
in  Eq,  F  is  given  by 


F  =  U2  *  U1  *  |IU22  -  2  U12  +  uni  <27> 

where 

U.  =  «  u.  »  (28a) 

and 

U..  s  «u.u.» -«u.»«u.»  :  i,j  =  1,2.  (28b) 

ij  i  J  i  J  J 

The  integrals  in  Uj  and  are  easily  performed  to  give 
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NaGE' 


(29) 


U 


1 


and 


UH  =  NG2E,2A2a2  (30) 

where  the  dimersionless  anisotropy  parameter  A  is  defined  by 

A2  =  2[(aj  -  a2)2  +  (o2  -  a3)2 +(Qj  -  q3)2]/(45  a2)  (31) 

and  a  =  (a.  +  +  a,)/3  is  the  usual  linear  molecular  polariza¬ 

bility.  1  C  5 

To  perform  the  integrals  in  U2  and  U  12  we  require  the 
positional  average  of  .  If  the  dimersionless  two-particle  cor¬ 
relation  junction  defineaoy 


p(rl2)  s  V  dr3  ...  dr^  w(jr  f ) 


(32) 


approaches  a  constant  value  when  rj?  =  |£l  “  £2  |  is  still  small  com¬ 

pared  with  the  ize  of  the  liquid  sample  under  consideration,  then  it  is 
well  known  that  (neglecting  terms  of  order  N"*  smaller)  the  required 
positional  average  gives  * ' 


>  -  "O'  <33> 

where 

I  =  L  —  4ir/3  ,  (34) 

1 7 

independently  of  the  form  of  p(r).  p  is  the  number  density  of  mole¬ 
cules  N/V,  L  is  the  depolarizing  factor  along  the  principal  axis  of  the 
ellipsoidal  sample  which  is  parallel  to  the  external  field.  Therefore 

U2  =  Nq2  E’2  pl/2  (35) 

and 

U12  =  Na3G  E'3  pi  A2  .  (36) 
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To  perform  the  integrals  required  in  U^Z*  t^ie  exact  f°rm  of 
the  two,  three,  and  four  particle  correlation  functions  is  required. 
For  purposes  of  calculation,  we  have  assumed  that  p(r)  is  a  simple 
step  function  which  properly  normalized  is 


where 


p(r)  =0,  r  <  d 

=  (1  -  t/V)"1,  r  >  d 


(37a) 


t  h  4ird3/3  . 


(37b) 


We  have  used  the  usual  "Kirkwood  superposition  approximation"  for 
the  three*  and  four-particle  correlation  functions,  expressing  them  in 
terms  of  p(r).  The  details  of  the  calculation  of  U22  with  these  approx¬ 
imations  is  given  in  Appendix  IV;  the  result  is 

U22/(NE'4  p2  a4)  =  (i  2  +  2tt2/3)  A2  +  2tt2/3 

+  tt2(32  +  136  A2  +  167  A4)/(45  p  r)  (38) 

+  0.961  p  t  (1  +  K4/4) 

where  K.  is  a  small  contribution  (~  10  *)  from  an  integral  that  we 
have  not  been  able  to  evaluate  exactly  and  which  we  may  neglect  without 
affecting  the  probable  accuracy  of  the  over -all  results. 

5.  The  Best  Effective  Field- Linear  and  Nonlinear 

Index  Formulae 


Having  explicit  expressions  for  all  the  terms  in  (27)  for 
F,  we  will  adjust  the  magnitude  E'  of  the  effective  local  field  to  obtain 
as  low  a  value  as  possible  for  the  upper  bound  F  of  the  free  energy. 
Since  the  terms  of  order  E^  are  supposed  to  be  orders  of  magnitude 
smaller  than  those  of  order  E&  ,  it  is  easily  seen  that  the  best  E'  is 
that  which  minimizes  the  terms  of  Eq.  Any  small  change  of  E'  from 
such  a  value  raises  these  large  terms  more  than  it  can  depress  the 
smaller  terms  of  order  E£.  That  is,  we  need  only  minimize  (U2  -Uj) 
with  respect  to  E'  and  find 

E'  =  Eo/(l  +  pal )  (39) 
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is  the  best  choice.  Using  (39)  in  ( 27 )  gives  immediately  the  approxi 
mate  b  and  c  coefficients  of  (7),  from  which,  with  the  aid  of  (7) 
through  (12),  the  following  relations  result  for  the  linear  dielectric 
constant  c  and  the  nonlinear  index  n^: 


^  1  +  Sirpa/3 

•  2  T~--firpa7T 


and 


n. 


2 

irpa  0 


___  \aZ  +  —  Pa - (32  +  136  A2  +  167  A4) 

n(l  -  4irpa/3)  l  T 

+  +  3.  843  a2  p3t]  . 


(41) 


We  reiterate  that  the  estimate  (41)  of  the  nonlinear  index  does 
not  contain  any  effects  of  macroscopic  density  change  or  current  in 
the  fluids,  but  presumably  does  represent  the  steady  state  effects  of 
all  microscopic  molecular  reorientation  and  redistribution. 


6. 


Discussion  and  Examination  of  Accuracy 


inequality  (40)  is  equivalent  to  i 
otti  function  (c  -  1 ) p ”  («  +  3)"1 


the  statement  that 


The 

the  Clausius -Mossotti  function  (c  -  l)p_1  («  +  3)"1  is  always  larger 
than  4ttq/3  fox:  a  classical  fluid  of  molecules  of  fixed  linear  polariza¬ 
bility,  regardless  of  the  intermolecular  forces,  provided  that  the  two- 
particle  correlation  function  depends  only  on  the  intermolecular  spacing 
and  becomes  constant  for  large  spacing.  The  fact  that  many  fluids  have 
been  found  which  slightly  disobey  (40)  demonstrates  that  at  least  one 
of  these  assumptions  is  not  entirely  valid.  However,  the  observed 
deviation  of  the  Clausius -Mossotti  function  from  4tt<i/3  is  rarely  more 
than  a  few  percent  and  this  suggests  that  our  nonlinear  term  (41)  might 
be  accurate  to  within  10  to  50>  ,  depending  on  the  molecule.  However, 
unlike  «,  is  very  sensitive  to  the  form  of  the  two-,  three-,  and 
four-particle  correlation  functions,  and  the  errors  that  resu.V  from 
approximating  the  latter  with  the  Kirkwood  approximation  are  completely 
unknown  at  present.  Nevertheless,  it  is  reassuring  that  certain  features 
of  (41)  are  familiar.  The  first  term  in  the  brackets  on  the  RHS  is 
exactly  the  formula  of  Debye1  .  for  the  dc  nonlinear  dielectric  constant 
(except  we  are  extending  it  to  optical  frequencies).  This  term  dominates 
all  others  at  low  densities.  .The  other  terms  are  all  proportional  to 
higher  powers  of  the  density  than  this  term  and  represent  high  density 
corrections  which  will  be  important  in  liquids. 
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Several  terms  in  (41)  require  knowledge  of  the  molecular 
volume  t/4.  This  parameter  also  occurs  in  the  theories  of  the  fluid 
viscosity,  the  heat  conductivity,  and  the  virial  coefficients.  An  ex¬ 
amination  of  data  on  representative  fluids  in  each  of  these  categories 
gives  values  for  t  which  vary  typically  by  10%  in  a  given  fluid. 

Data  are  not  available  for  all  of  the  fluids  of  interest  to  us,  but  in 
those  where  it  is  available,  t  is  always  within  ±  10%  of  2.  5/  p. 
Therefore,  we  have  used  this  value  for  t  for  all  liquids  in  numerical 
calculations  of  the  indices.  The  nonlinear  indices  nz  predicted  by 
(41)  are  listed  in  Table  I  for  those  liquids  whose  nonlinear  effects 
have  been  commonly  studied.  The  table  also  includes  several  sym¬ 
metric-molecule  liquids  whose  nonlinear  indices  are  predicted  to  be 
comparable  with  the  most  nonlinear  asymmetric-molecule  liquids. 
However,  the  molecules  in  these  liquids  are  likely  to  be  so  large 
that  our  model  of  hard  spheres  interacting  as  if  point  dipoles  were 
at  their  centers  breaks  down.  The  values  of  nz  for  these  molecules 
are  therefore  intended  to  be  suggestive  of  possible  effects. 

With  the  improvement  in  cur  approximations  when  applied  to 
the  smaller  molecules  and  from  comparing  our  results  with  avail¬ 
able  data,  we  expect  that  our  tabulated  values  for  n z  (for  these 
smaller  molecules) should  have  a  relative  accuracy  of  perhaps  ±  20% 
or  better.  The  absolute  accuracy  is  difficult  to  judge  because  of  the 
uncertainties  in  the  Kirkwood  approximation. 


Even  considering  the  large  “possible  errors  in  (41),  one  is  led  to 
conclude  that  reorientation  effects  on  the  index  are  not  much  more 
important  than  redistribution  effects,  and  often  contribute  negligibly 
to  the  change  in  index. 

In  addition,  the  recent  experimental  results  of  Mayer  for  two 
thirds  of  the  difference  between  the  nonlinear  indices  parallel  and 
perpendicular  to  the  applied  field*-*  are  given  in  the  table.  If  mole¬ 
cular  reorientation  were  the  dominant  mechanism  and  the  Debye 
theory*®  correct,  this  number  would  equal  n^. 

Also  included  in  Table  I  are  the  recent  experimental  results 
of  Wang  and  Racette^  for  a  coefficient  Cj£21  which  expresses  the 
nonlinear  tendency  of  the  liquid  to  rotate  the  axes  of  the  polarization 
ellipse  of  the  incident  light.  If  molecular  orientation  were  dominant 
and  the  Debye  theory*®  correct,  the  numbers  in  this  column  would 
also  equal  n^. 

Unfortunately,  our  theory  of  the  nonlinear  index  is  only  for 
the  change  in  index  for  light  polarized  parallel  to  a  strong  plane 
polarized  beam.  When  the  theory  is  extended  to  two  beams  of  different 
polarization,  the  ideas  about  molecular  redistribution  can  be  checked 


Included  in  Table  I,  along  with  our  estimates  of  n^  and  Pc, 
is  the  nonlinear  index  n?  which  would  result  from  the  low  density 
theory  of  Debye  which  neglects  the  effects  of  molecular  redistri 
bution.  The  Debye  index  n?  is  always  smaller  than  that  of  our  (41] 
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TABLE  I 


Some  Nonlinear  Optical  Coefficients  for  21  Liquids* 


Li  quid 

0  ,n12 
n^  x  10 

calc,  from 
Ref.  18 

i  nl  2 
n^  x  10 

calc,  from 
eq.  (47) 

T  '  "2i> 

measured , 

Ref.  15 

10 1  ^  x  16trc  1 22  j/n 

measured 

Ref.  19 

Pc(kW)  calc, 
for  6943  ft 

1. 

Lead  tetraisopropyla 

0 

111 

1.51 

1 

Tin  tetra  2-metholbutyla 

0 

109 

1. 54 

1. 

Tin  tetrapental* 

0 

108 

1. 55 

4. 

CS2 

22 

96 

13.6 

10.5 

1. 74 

5. 

Tin  tatralbutyla 

0 

89 

1.88 

6. 

Si  Br* 

4 

0 

82 

2.  05 

7. 

Lead  tetraethyl* 

0 

81 

2.  07 

8. 

Benzoylrhloride' 

73b 

7.9 

2.  31 

9. 

Nitrobenzene 

8.  8 

69 

9.4 

6.  8 

2.  43 

10. 

Bromobenzene 

6.  4C 

68 

2.47 

11. 

Acetophenone 

56d 

3.  10 

12. 

Chlorbenzene 

6.  3 

54 

3.  1  1 

13. 

Toluene 

5.  6 

45 

3.  0 

2.  7 

3.77 

14. 

Benzene 

4.  4 

37 

1.3 

2.2 

4.5 

15. 

c  ci4 

0 

24 

1.7 

7.  1 

lb. 

Chloroform 

1.  3 

19 

0.  54 

9.0 

17. 

Water 

llb 

0.  1 

15.4 

18. 

Acetone 

0.  6 

6.  7 

0.  2 

25.5 

19. 

Acetic  acid 

5.  7b 

0.  4 

30 

20. 

Liquid  CH^ 

0 

3.  0 

56 

21. 

Liquid 

0.  010 

0.  024e 

7000 

♦  The  values  of  n£  are  calculated  from  the  standard  relation  of  (18)  for  the  effect  of  reorientation  of  asymmetric 
molecules.  The  n2  are  calculated  from  (47)  and  include  redistribution  as  well  as  reorientation  effects.  The 

2/ 3(n2  -  n2x)  are  t*1*  values  of  optically  induced  birefrigence  measured  by  Mayer*5  and  would  equal  n2  if  the 
reorientation  theory18  were  valid.  The  l6wCi22i/n  were  measured  by  Wang  and  Racette1^  and  would  equa( 
n2  if  the  reorientation  theory  were  valid.  The  critical  powers  for  self-trapping  are  computed  using  our  calcu¬ 
lated  values  for  n^  in  (12).  The  liquids  are  given  in  the  nominal  order  of  the  Pc.  The  accuracy  of  the  given 
results  is  discussed  in  the  text  and  is  less  than  indicated  by  the  figures  given  here.  The  polarizability  data 
used  in  the  computations  are  from  Ref.  20. 

Purely  symmetric  molecule. 

Estimated  by  adding  value  from  column  3  to  n,  calculated  with  A  =  0. 

""Estimated  by  scaling  Table  142072  in  Ref.  20  to  Table 

142071. 

Estimated  by  assuming  molecule  is  symmetric  (A  =  0] 

♦ 

This  value  is  less  accurate  than  others  in  column  because  exhibits  quantum  effects  which  we  have  ignored. 
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by  direct  comparison  with  the  results  of  Mayer  and  of  Wang  and 
Racette.  * '  This  extension  is  planned. 

7.  Conclusions 


The  results  of  (41)  as  expressed  in  Table  I  indicate 
that  no  room  temperature  liquid  can  be  expected  to  support  a  coher¬ 
ent  beam  of  as  much  as  a  megawatt  for  many  diameters  of  path 
length  without  significant  beam  deterioration  occurring.  Even  liquid 
hydrogen,  which  has  a  much  lower  expected  nonlinear  index  than 
normal  liquids,  would  not  be  expected  to  support  beams  of  much  over 
107  Wr  without  self-focusing  effects.  However,  we  note  the  following 
points  which  may  bear  on,  or  modify,  the  predicted  critical  powers 
of  Table  I  and  hence  may  modify  the  conclusions. 

First,  the  first  quantitative  measurements  of  power  thres¬ 
holds  for  self- trapping  in  benzene,  nitrobenzene,  and  toluene',  and 
in  CS?  (Ref.  10)  have  given  thresholds  an  order  of  magnitude  greater 
than  those  predicted  in  Table  I.  However,  Garmire,  Chiao,  and 
Townes  (private  communication)  have  recently  seen  "small  scale" 
trapping  which  sets  in  at  an  as  yet  undetermined  power  level  some¬ 
what  lower  than  that  which  they  reported  in  Ref.  (10).  Therefore, 
the  size  of  the  discrepancy,  if  any,  between  our  predicted  thresholds 
and  those  observed  is  still  uncertain. 

Secondly,  as  Bloembergen  and  Lallemand  have  pointed  out, 
the  time  delay  for  some  molecules  to  reorient  themselves  in  a 
suddenly  applied  large  electric  field  may  not  be  small  compared  with 
the  typical  duration  of  a  laser  pulse.  This  suggests  that  if  the  powers 
used  in  the  previous  measurements  of  thresholds  were  kept  on  longer, 
the  thresholds  might  have  been  lower.  That  is,  longer  laser  pulses 
may  be  required  if  the  "steady  state"  nonlinear  index  we  have  calcu¬ 
lated  is  to  be  achieved.  Crude  estimates  of  this  relaxation  time  made 
on  the  basis  of  the  classical  Debye  theory  indicate  that  liquids  having 
a  viscosity  greater  than  or  of  the  order  of  one  poise  might  encounter 
this  raising  of  threshold.  The  liquids  of  Table  I  are  less  viscous  than 
this  by  an  order  of  magnitude  or  more. 

Third,  stimulated  Brillouin  scattering  may  have  a  threshold 
lower  than  that  for  self- trapping  and  complicate  matters  in  a  way 
which  spoils  further  the  coherence  of  any  Raman  laser  action. 

Also,  the  question  of  how  the  turbulence  and  vibrational  waves 
which  exist  in  liq-  'ds  under  the  normal  conditions  in  which  experi¬ 
ments  are  perfom.  d,  affects  the  threshold  for  self-focusing  has 
never  been  answered.  Small  bubbles  and  impurities  may  also  have 
some  effect.  All  such  effects  would  probably  raise  the  threshold  for 
producing  observable  self-focusi  g  above  the  apparently  low  predic¬ 
tions  in  Table  I  of  (41). 
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IV. 


SUMMARY 


The  self-focusing  or  self-trapping  of  the  pump  light  or  the 
Raman  scattered  light  in  a  Raman  laser  tends  to  destroy  the  coher¬ 
ence  of  the  Raman  output.  Therefore,  self-focusing  effects  must 
be  understood  and  probably  prevented  if  a  high  coherence,  high 
energy  Raman  output  is  to  be  obtained.  Our  experimental  studies 
of  stimulated  Raman  scattering  and  the  concomitant  beam  trapping 
have  demonstrated  the  following  points. 

Although  the  time-integrated  intensities  of  observed  fila¬ 
ments  formed  in  the  laser  beam  by  a  Ramar.-active  liquid  are  not 
nearly  high  enough  to  explain  the  observed  anomalously  low  SRS 
thresholds,  the  instantaneous  intensities  we  have  observed  by 
time- resolved  photographic  techniques  are  indeed  high  enough.  A 
given  ruby  and  Raman  filament  never  seems  to  persist  for  more 
than  a  small  fraction  of  the  laser  pulse  duration;  different  filaments, 
however,  occur  at  different  times  during  the  pulse.  The  "gain 
anomaly"  is  therefore  resolved,  at  least  for  the  present,  in  the 
commonly  studied  liquids  —  CS^  ,  benzene,  nitrobenzene,  and 
toluene. 

If,  as  is  commonly  supposed,  a  molecular  reorientation 
were  mainly  responsible  for  the  nonlinear  index  and  for  self- trapping 
in  the  Raman-active  liquids,  a  simple  theory  we  developed  suggests 
that  a  circularly  polarized  plane  input  wave  should  exhibit  a  thres¬ 
hold  for  self-trapping  four  times  that  of  a  linearly  polarized  wave. 
Our  experimental  efforts  to  check  this  ratio  showed  that,  in  fact, 
a  stable  mode  of  propagation  does  not  exist  for  any  state  of  polari¬ 
zation  except  linear  in  CS£  ,  nitrobenzene,  and  presumably  many 
other  liquids  as  well.  As  soon  as  self- trapping  in  a  circularly 
polarized  beam  was  observed,  the  trapped  light  was  seen  to  be  de¬ 
polarized.  We  have  explained  this  behavior  qualitatively,  and 
expect  it  to  occur  for  any  known  physical  mechanism  which  can 
cause  a  nonlinear  index,  except  electrostriction.  We  therefore 
conclude  that  electrostriction,  often  proposed  as  a  mechanism  for 
self- trapping,  is  unlikely  to  be  a  major  factor  in  self- trapping  in 
many  liquids.  However,  we  were  not  able  to  conclude  anything 
about  the  reorientation  hypothesis  from  our  observations. 

Our  theoretical  studies  of  beam  trapping  and  its  relation  to 
SRS  have  been  mainly  devoted  to  understanding  how  self- trapping 
could  be  eliminated  or  at  least  made  tolerable.  To  this  end  it 
was  important  to  ascertain  the  relative  importance  of  various  physi¬ 
cal  mechanisms  in  creating  a  nonlinear  index  in  various  situations. 
We  devised  a  3elf-trapping  threshold  ratio  test  and  computed  the 
expected  ratios  for  different  mechanisms.  When  observed  insta¬ 
bility  in  the  propagation  of  nonlinearly  polarized  light  made  this  test 
infeasible,  we  could  only  infer  that  electrostriction  was  not  a 
dominant  mechanism. 
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We  therefore  made  a  statistical  mechanical  calculation  of 
the  nonlinear  index  which  assumed  that  no  macroscopic  density- 
changes  occurred,  such  as  might  occur  from  electros tricti on  or 
heating.  It  was  hoped  that  molecules  with  no  asymmetry  of  polari¬ 
zability  would  shov  only  the  very  small  electronic  ^onlinearities 
in  their  dielectric  behavior.  This  hope  proved  in  vain,  because  the 
calculations  showed  that  molecular  redistribution  in  space  would 
cause  nonlinearities  comparable  with  those  caused  by  reorientation 
of  very  asymmetric  molecules.  We  conclude,  therefore,  that  no 
room  temperature  liquid  will  ever  be  likely  to  sustain  a  high  energy 
Raman  output  beam  without  its  nonlinear  dielectric  properties 
significantly  distorting  the  beam.  We  look  to  gases  as  the  most 
promising  avenue  for  further  investigatici*  of  the  high  energy,  high 
coherence  potential  of  Raman  lasers. 
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Longitudinal  Mode  Control  in  Giant  Pulse  Lasers 

F.  J.  McCLUNG  and  D.  WEINER 


Abstract—  W*  have  introduced  longitudinal  mode  selection 
into  a  (iant  pulse  laser  to  obtain  single  mode  output  from 
the  laser.  Some  advantages  of  achieving  longitudinal  mode 
control  in  a  giant  pulse  laser  are  noted.  The  methods  of  mode 
control  used  are  described.  These  include  orienting  various  re¬ 
flecting  surfaces  in  the  laser  cavity,  cooling  the  ruby  laser  crystal, 
and  introduction  of  a  saturable  dye  in  the  cavity.  Methods  of 
measuring  the  mode  structure  are  given.  Results  of  these  meas¬ 
urements  with  varying  degrees  of  mode  control  are  described  in 
detail.  Under  some  conditions  it  is  possible  to  obtain  essentially 
single  mode  behavior  from  the  giant  pulse  laser. 

Introduction 

The  FREQUENCY  SPREAD  of  the  output  of  a 
rotating  prism  giant  pulse  ruby  laser  (GPL)  is 
typically  about  one  Angstrom  at  room  tempera¬ 
ture.  Closer  examination  shows  this  output  to  consist 
of  many  spectrally  sharp  components.  These  components 
are  spaced,  as  may  be  expected,  by  frequencies  equal  to 
the  reciprocal  of  twice  the  transit  time  between  various 
pairs  of  reflecting  surfaces  in  the  laser  cavity. 

For  many  reasons  it  would  be  advantageous  to  reduce 
the  spread  in  frequencies  as  much  as  possible  so  that  ul¬ 
timately  only  a  single  longitudinal  mode  of  the  laser 
cavity  is  operating.  This  would  make  it  possible  to  ob- 
M'vve  -mall  frequency  shifts  in  experiments,  such  as 
Hrillouin  scattering  or  Rayleigh  scattering.  If  the  total 
output  can  he  held  essentially  constant  (as  has  been 
found  possible) .  one  would  have  a  higher  spectral  bright - 
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ness.  This  in  turn  would  imply  enhanced  nonlinear  ef¬ 
fects  in  systems  whose  natural  frequency  spread  is 
smaller  than  that  of  the  GPL  as  would  be  the  case  for 
stimulated  Raman  scattering  (SRS)  in  gases  or  for  non¬ 
linear  optical  effects  in  low-density  plasmas.  Also,  the 
comparison  of  theory  and  experiment  is  much  easier  for 
a  single  mode  laser  pump  for  such  processes  as  SRS. 
Here,  for  example,  multimode  theory  [1]  is  qualitative 
in  predicting  the  complex  characteristics  of  anti-Stokes 
radiation  for  a  multimode  pump.  Finally,  it  is  sometimes 
|K>ssible  to  observe  deep  ( ~50  per  cent)  temporal  modu¬ 
lation  of  the  GPL  pulse  corresponding  to  mode  spacings 
of  the  cavity.  Thus,  mode  control  techniques  offer  in¬ 
sights  into  methods  of  pulse  shaping  at  these  high  fre¬ 
quencies. 

Wc  have  achieved  single  mode  behavior  by  a  combina¬ 
tion  of  techniques.  These  include  carefully  orienting  vari¬ 
ous  reflecting  surfaces  in  the  laser  cavity,  cooling  of  the 
ruby  laser  crystal,  and  introduction  of  a  saturable  dye 
in  the  cavity.  We  shall  describe  these  techniques  in  de¬ 
tail  and  also  discuss  methods  of  measuring  the  mode 
structure  in  what  follows. 

Techniques  of  Mode  Control 

Figure  1  shows  a  schematic  diagram  of  the  apparatus 
used  for  longitudinal  mode  selection.  The  choice  and 
arrangement  of  the  elements  are  by  no  means  unique; 
they  are  in  part  determined  by  such  considerations  as 
the  availability  of  optical  flats  of  suitable  thickness  and 
index  of  refraction,  the  availability  of  suitable  optical 
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!'V  1.  An  experimental  arrangement  for  longitudinal  mode  con¬ 
trol  of  a  giant  ptiW  laser  showing  gas  laser  (GL),  prism  (P), 
heiiin  splitter  (B5),  reflection  mode  selector  (RMS),  transmis¬ 
sion  inode  selector  (TMS),  ruby  (R),  flashlamp  housing  (FH), 
crvptocynnine  mode  selector  (CMS),  rotating  prism  (RP), 
thermistor  (Tl),  thermocouple  (T2),  resistance  bridge  (RB), 
voltmeter  (VM),  valve  (V),  liquid  nitrogen  (LN),  compressed 
nitrogen  (CN),  and  Dewar  (D). 

coatings  and  the  constraints  imposed  by  the  dimensions  of 
the  existing  GPL.  We  make  no  effort  here  to  determine  the 
influence  on  our  results  of  deviation  from  flatness  or 
parallelism  in  the  various  optical  elements.  The  ruby  is 
usually  poor  (2a  at  the  best  for  our  ruby),  and  our  opti¬ 
cal  flats  arc  finished  to  approximately  A/10  and  parallel 
in  transmission  to  A/10. 

In  the  future  this  inode-selected  laser  may  well  be 
made  simpler.  However,  the  present  apparatus  is  suffi¬ 
cient  to  produce  essentially  single  mode  output  with  a 
power  approximately  that  of  a  nonmode-sclcctcd  GPL. 
It  also  serves  well  to  illustrate  the  various  methods  for 
mode  selection. 

Reflection  Mode  Selection 

Without  longitudinal  mode  selection  a  typical  GPL 
has  a  linewidth  of  ~1  —  cur1  which  correcponds  to 
— 150  modes  for  a  50-cm  laser  cavity.  By  replacing  the 
front  reflector  of  the  GPL  with  a  suitable  optical  flat  we 
may  modulate  the  mode  structure  so  as  to  greatly  reduce 
the  number  of  oscillating  modes  as  first  reported  by 
Stoicheff  [2]. 

Our  reflection  mode  selector  (RMS  in  Fig.  1)  is  a 
— 2-mm-thick  interferometrieally  flat  and  parallel  sap¬ 
phire  plate  (n  =  1.75).  The  spacing  between  modes  of 
this  flat  is  1.4  cm-1.  Hence,  if  one  reflection  maximum  is 
centered  on  the  peak  of  the  R,  line  in  ruby,  the  gain  at 
the  adjacent  reflection  maxima  should  be  insufficient  to 
produce  oscillation.  In  this  way  the  number  of  possible 
cavity  modes  is  reduced  from  —  150  to  —20. 

The  reflection  mode  selector  is  aligned  parallel  to  the 
uncoated  ruby  (R)  and  rotating  prism  Q  switch  (RP) 
with  the  aid  of  a  gas  laser  alignment  light  source  (GL) 
whose  beam  enters  the  laser  cavity  via  a  prism  (P)  and 
a  beam  splitter  (BS).  As  far  as  mode  selection  is  con¬ 
cerned,  alignment  of  the  surfaces  by  superposition,  in  the 
far  field,  of  reflections  from  the  various  surfaces  is  al¬ 
most  as  good  as  alignment  by  observing  interference 
effects  in  reflected  light  from  the  gas  laser.  However,  the 
latter  method  seems  to  give  a  higher  power  output. 


ima  of  the  reflection  mode  selector  by  cooling  the  ruby. 
In  our  apparatus  the  ruby  is  cooled  to  about  0°C.  This 
cooling  is  advantageous  in  other  ways.  It  produces  a 
decrease  in  the  width  of  the  Rx  line  which  further  reduces 
the  number  of  possible  modes.  It  also  increases  the  peak 
gain  for  the  R i  line.  This  is  useful  because  the  peak  re¬ 
flectivity  of  the  reflection  mode  selector,  whose  refractive 
index  was  as  high  as  was  readily  available,  was  only 
—25  per  cent,  without  the  added  reflectivity  from  the 
ruby. 

The  cooling  is  accomplished  by  passing  cold  dry  nitro¬ 
gen  gas  past  the  ruby  crystal.  Compressed  dry  nitrogen 
(CN)  is  passed  into  the  flash  lamp  housing  (EH)  and 
past  the  ruby.  Part  of  this  nitrogen  has  been  previously 
diverted  through  a  heat  exchanger  that  is  in  contact  with 
liquid  nitrogen  (LN)  in  a  Dewar  (D).  The  amount  of 
cooled  nitrogen  and  hence  the  temperature  of  the  gas 
mixture  reaching  the  ruby  is  controlled  by  a  valve  (V). 
The  inlet  temperature  of  the  gas  is  monitored  by  a 
copper-constantan  thermocouple  (T2)  whose  voltage  is 
read  on  a  microvoltmeter  (VM).  The  temperature  of  the 
ruby  is  monitored  by  a  thermistor  (Tl)  whose  tempera¬ 
ture  is  determined  by  the  current  imbalance  in  a  resist¬ 
ance  bridge  (RB). 

We  wish  to  temperature- tune  the  ruby  so  that  its 
fluorescence  peak  coincides  with  a  reflection  maximum 
of  the  sapphire-plate  mode  selector.  Near  room  tempera¬ 
ture,  we  can  accomplish  this  by  observing  the  mode 
structure  with  a  Fabry-Perot  interferometer  at  a  series 
of  ruby  equilibrium  temperatures  spaced  1°  or  2°  apart. 
For  two  of  these  temperatures  differing  by  about  10°, 
the  output  of  the  GPL,  measured  with  low  resolution, 
consists  of  two  modes  of  equal  intensity  spaced  1.4  cm-1 
apart.  The  desired  operating  temperature  is  halfway  be¬ 
tween  these  two  temperatures.  At  temperatures  close  to 
0°C,  it  is  better  to  determine  the  proper  temperaure  by 
high  resolution  observation  of  the  GPL  output. 

To  reduce  the  number  of  oscillating  modes  still  fur¬ 
ther,  one  may  introduce  other  sets  of  interferometrieally 
flat  surfaces  within  the  laser  cavity.  These  have  trans¬ 
mission  maxima  and  minima  whose  frequency  spacing 
is  determined  by  the  distance  between  the  surfaces  and 
by  the  refractive  index  of  the  intervening  material. 

The  ruby  crystal  provides  one  such  set  of  surfaces. 
Our  3-in. -long  (f-in.  diameter)  ruby  has  transmission 
maxima  every  1.2  Gc.  Thus  it  will  suppress  two  or  three 
out  of  every  four  consecutive  cavity  modes  that  are  sepa¬ 
rated  by  300  Me.  What  remains  then  is  to  suppress  all 
but  one  of  the  approximately  five  modes  of  the  ruby 
which  are  separated  by  1.200  Gc.  This  may  be  done,  in 
principle,  by  placing  the  ruby  at  such  a  distance  from 
another  reflecting  surface  so  as  to  form  a  suitable  cavity. 
This  may  well  have  been  done,  intentionally  or  other¬ 
wise,  in  other  experiments  where  mode  control  has  been 
achieved. 

The  total  effect  of  placing  sets  of  reflecting  surfaces 
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I *;i v:i  1  lt-1  in  oilier*  in  l hr  laser  cavity  i*.  t»l"  course.  very 
«*« >in| di'x.  Each  surface  form*  :i  resonator  with  any  other 
-url'ace.  thus  giving  a  lnultitude  of  interacting  resonant 
ravine*  of  different  frequencies.  The  net  olltjilll  spectrum 
of  the  sy*trm  will  also  depend  upon  the  reflectivities  of 
th<-  vnrioii*  *itrl'aee.»  ami  the  gain  ol'  any  material  in  the 
caviu  In  otiler  to  gain  more  intuition  concerning  the 
I'l'ohleiu.  Wagner  |3|  of  this  laboratory  performed  eom- 
I'liter  ealeiilation*  to  lind  the  n linin tin n  gain  necessary 
for  o*eillation  for  many  consecutive  modes  in  t  system 
eon-Ming  of  a  roof  prism,  a  rnhy.  a  glass  lelleel ion  mode 
*rlrelor,  and  a  l’i-per  rent  end  reflector.  These  ealrnla- 
tion*  si io\\*  that  the  position  of  the  rnhy  ami  the  reso¬ 
nance*  of  the  roof  prism  are  relatively  unimportant.  They 
al*o  show  1 1  wit.  the  re*onanees  of  the  rnhy  (perhaps  due 
to  it*  gainl  and  the  position  and  thickness  of  the  «ilnss 
rdleetion  mode  sector  are  important.  These  observa¬ 
tions  help  in  choosm*:  onr  initial  cavity  con  figuration. 
However,  in  onr  present  configuration  (Fig.  1)  we  find 
that,  the  position  of  the  rnhy  is  important  as  it  forms, 
with  the  reflection  mode  selector,  ail  effective  reflector  of 
-nffieient  reflectivity  to  produce  oscillation. 

Tnuisnilsxiiw  Moth  Selection 

Heeause  of  the  physical  limitations  in  onr  existing 
(11*1..  it  is  not  possible  to  position  the  rnhy  to  suppress 
all  hut  one  of  the  1.2-(5e  modes.  To  perform  this  task  we 
in-ert  an  optically  flat  and  parallel  glass  plate  that  is 
tilted  slightly  with  respect  to  the  laser  beam.  We  can 
understand  the  operation  of  this  mode  selector  by  coit- 
sidering  that  this  tilted  plate  has  no  reflection  lo*>  for 
frequencies  corresponding  to  its  Fnbrv-Perut  transmis¬ 
sion  maxima.  At  other  frequencies  the  reflections  from 
this  mode  selector  are  lost  from  the  cavity  and  thus  con¬ 
stitute  a  frequency-dependent  loss  mechanism  at  a  given 
angle.  The  fact  that  the  plate  is  tilted  means  that  there 
me  no  other  resonators  formed  with  other  reflectors  in 
tie  cavity.  The  frequency  is  tuned  by  varying  the  tilt 
of  tiie  flat. 

The  transmission  mode  selector  (T.MSi  used  is  a  1-cin- 
thiek  interferometrieally  flat  piece  of  flint  glass  whose 
modes  are  separated  by  —  10  Ge.  The  correct  angle  of 
tilt,  is  determined  empirically,  i.c..  by  seeing  which  of 
several  angles  gives  the  least  number  of  1.2-Ge  modes 
and  mode  interaction  due  to  the  1-rm  flat  Tl.e  relative 
angle  of  tilt  is  determined  by  observing  the  deviation 
angle  of  the  gas  laser  beam  upon  reflection  from  the  flat. 
The  use  of  this  mode  selector  also  has  the  advantage  of 
'Oppressing  weak  side  inodes  of  the  reflection  mode  sc- 
lector  which  otherwise  could  have  -ppeaml  at  high 
pumping  levels. 

Dili  Motlc  Selector 

At  the  high  power  levels  desired  f  —  1 0  MW  I,  the  fore¬ 
going  techniques  may  prove  insufficient  to  produce  re¬ 
liable  *ingle  mode  behavior.  Hence,  one  may  add  a  cell 
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containing  a  eryptoeyanine  dye  soluiion  in  the  foregoing 
configuration  as  had  been  shown  useful  in  mode  control 
j 4 1 .  The  eryptoeyanine  mode  selector  (C'.MSi  we  u*e  is 
formed  by  a  —  5  X  lb  molar  solution  of  eryptoeyanine 
in  nietlianol  placed  in  a  cell  with  interferometrieally  flat 
windows  2  mm  apart.  Its  absorption  at  G943  A  i*  —23 
per  cent.  The  cell  is  wedged  and  tilted  slightly  with  rc- 
*pcct  to  the  laser  beam  to  reduce  trnnsmi-sion  mode 
selection  effects. 

Measurement  of  Mode  Sthictirf. 

Our  primary  method  of  observing  the  mode  structure 
of  onr  GPI.  is  with  Fabry -Perot  interferometer*.  Fnr 
low-resolution  studies  we  use  »  2.4-nun  separation  be¬ 
tween  interferometer  plates.  For  high -resolution  work 
the  spacing  is  20  nun  between  plates  which  gives  an  in- 
terordcr  separation  of  7.3  Ge.  Each  plate  is  flat  and 
parallel  to  1/80  wavelength,  is  coated  with  silver  to  a 
transmission  of  —  1  per  cent,  and  is  —2  cm  in  diameter. 
The  instrument  is  used  with  a  camera  of  1-m  foeal 
length.  The  resolution  with  the  2-cm  spacer  is  —200  Me. 

Some  information  about  the  mode  structure  can  be 
obtained  from  the  temporal  modulation  of  the  GPL  out¬ 
put.  This  is  observed  with  a  biplanar  photodiode  and  a 
Tektronix  519  oscilloscope  of  combined  resolution  —  1 
Ge.  Under  certain  circumstances  we  observe  a  —50  per 
cent  modulation  at  the  fundamental  cavity  frequency. 

Stoieheff  [2]  has  also  shown  that,  the  spectral  width 
of  the  anti-Stokes  radiation  in  stimulated  Raman  scat¬ 
tering  (HRS)  is  very  sensitive  to  mode  structure  and  that 
very  narrow  lines  are  to  he  expected  with  single  longi¬ 
tudinal  mode  behavior.  Such  behavior  has  been  observed 
in  our  laboratory  and  is  shown  in  Figs.  2  and  3.  This  is 
perhaps  our  most  sensitive  test  for  single  mode  behavior. 

It  has  lieen  suggested  f  1  ] ,  [51  that  the  conversion  effi¬ 
ciency  of  GPL  to  Stokes  radiation  in  SRS  is  also  sensi¬ 
tive  to  the  number  of  inodes  of  the  GPL  and  that  the 
conversion  efficiency  for  a  single  inode  pump  should  he 
about  an  order  of  magnitude  less.  Measurements  with 
our  mode-controlled  GPL,  whose  output  consists  of  two 
TEMmoh  inodes  spaced  by  800  Me  (one  of  which  was 
—  4  times  stronger  than  the  other)  showed  that  this  was 
not  the  ease.  On  the  contrary,  the  conversion  efficiency 
at  a  given  power  density  was  precisely  independent  of 
the  number  of  modes  operating. 

The  transverse-mode  stmeture  was  evaluated  by  ob¬ 
serving  the  far-field  pattern  of  the  laser  with  a  1-m-focal 
length  camera.  The  pattern  observed  with  transverse¬ 
mode  selection  consisted  of  one  lobe  (TEMoo)  or  two 
lobes  (TEMoi).  The  size  of  each  lobe  corresponded  to 
the  diffraction  limit  for  the  laser  beam  diameter.  These 
patterns  looked  exactly  like  those  observed  by  Evtuliov 
and  Neeland  [5]  for  single-transversc-mode  operation  in 
normal  ruby  lasers.  The  power  output  was  — 2MW.  Fur¬ 
ther  details  of  transverse-mode  selection  will  be  pub¬ 
lished  elsewhere. 
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Fife.  2.  First  anti-Stokes  stimulated  Hainan  scattering  in  benzene 
allowing  a  pliotogia|i|i  of  emission  angle  vs,  wavelength  without 
mode  selection.  I  lie  neon  reference  ones  are  separated  by  —26 
A.  This  distance  corresponds  to  •*-1.5°  on  the  vertical  scale.  Neu¬ 
tral  Density  (X.D.)  1  tilt  r  covers  half  the  slit.  Copper  shields 
were  placed  over  the  <•  ntcr  of  the  slit  where  the  laser  was 
focused. 


Fig.  3.  First  anti-Stokes  stimulated  Raman  scattering  in  benzene 
showing  a  photograph  of  emission  angle  vs.  wavelength  with 
mode  selection.  The  neon  reference  lines  are  separated  by  —26 
A.  This  distance  corresponds  to  —1.5°  on  the  vertical  scale. 
An  N.D.  1  filter  covers  half  the  -lit.  Copper  shields  were  placed 
over  the  center  of  the  slit  where  the  laser  was  focused. 

Results 

The  results-  of  the  forepoinjr  modi*  selection  procedure 
can  be  seen  in  the  various  Fabry -Perot  photographs  of 
the  laser  output  shown  in  Fig-.  4  8.  For  Fig.  4,  the  2.4- 
tmn-spaccd  ctalon  is  used,  and  for  the  others,  the  20-mm- 
spaced  etalon  is  used.  In  each  picture  a  Neutral  Density 
0.6  filter  covers  about  half  the  pattern. 

Figure  4  shows  the  result  of  using  the  temperature 
tuning  and  the  reflection  mode  selector  only.  The  photo¬ 
graph  shows  an  output  consisting  of  two  components 
separated  by  —3  Gc.  The  spectral  width  here  is  limited 


by  the  1-Gc  resolution  of  the  interferometer.  When  the 
temperature  is  raised  by  — 5°C,  we  observe  two  com¬ 
ponents  of  about  equal  intensity  separated  by  40  Gc. 
This  separation  corresponds  to  adjacent  modes  of  the 
2-mm-thick  reflection  mode  selector.  The  laser  output 
was  12  MW. 

Figure  5  shows  the  result  of  adding  our  transmission 
mode  selector  to  the  foregoing  configuration  and  adjust¬ 
ing  its  tilt  for  minimum  laser  linewidth.  The  laser  output 
is  —5  MW.  Careful  examination  of  the  photographic 
plate  shows  that  <-80  per  cent  of  the  energy  is  in  a  single 
component  —450  Me  wide,  which  probably  corresponds 
to  two  components  of  roughly  equal  intensity  spaced  by 
—300  Me.  The  rest  of  the  energy  is  in  four  weak  com¬ 
ponents  that  are  spaced  by  —1.2  Me.  We  note  that,  be¬ 
cause  of  halation  and  lack  of  dynamic  range,  such  an 
output  would  appear  to  be  a  single  longitudinal  mode  on 
Polaroid  film. 

Figure  6  shows  the  effect  of  adding  the  cryptocyanine 
dye  cell  to  the  foregoing  configuration  while  keeping  the 
pumping  level  constant.  The  output  is  —1  MW  and 
appears  as  a  single  component  —300  Me  wide,  which 
corresponds  to  a  single  longitudinal  mode. 

We  have  found  the  cryptocyanine  dye  insufficient  by 
itself  to  suppress  components  40  Gc  apart  when  the  laser 
output  is  — 10  MW.  At  these  output  levels  the  dye,  usu¬ 
ally  but  not  reliably,  suppresses  some  of  the  more  closely 
spaced  (—1  Gc)  modes  when  used  with  the  transmission 
mode  selector.  This  effect  is  shown  in  Figs.  7  and  8,  where 
the  output  is  maintained  at  —8  MW  by  pumping  the 
laser  harder  when  the  dye  is  used.  The  configuration  is 
somewhat  different  from  that  for  Figs.  5  and  6  in  that 
the  distance  between  the  reflection  mode  selector  and  the 
front  of  he  ruby  is  changed  from  12  cm  to  18  cm  and 
the  total  cavity  length  is  chang'd  from  50  cm  to  70  cm. 
Figure  7  shows  the  laser  output  without  the  dye.  There 
are  one  strong  and  four  weak  components,  each  —300 
Me  wide  and  spaced  by  —800  Me  from  the  adjacent 
components.  Figure  8  shows  the  effect  of  adding  the  dye. 
There  are  only  two  components  —300  Me  wide  and  of 
comparable  intensity  spaced  by  —800  Me.  The  800-Mc 
splitting  corresponds  to  the  new  ruby-to-reflector  dis¬ 
tance. 

Finally,  we  note  that  another  3-in.  X  f-in.  ruby  gives 
essentially  the  same  results  as  the  one  used  previously. 

Summary 

It  is  possible  to  mode-select  the  rather  broadband  out¬ 
put  of  the  ruby  giant  pulse  laser  without  serious  loss  of 
power.  The  spectral  purity  of  the  output  will  depend  to 
some  extent  upon  the  power  output  desired  from  the 
laser.  In  general,  it  appears  possible  to  obtain  as  much 
as  80  per  cent  of  the  energy  in  a  single  longitudinal  mode 
of  the  laser  without  sacrificing  more  than  30  per  cent  in 
the  output  power.  A  comparison  of  the  observed  line- 


I' in.  I.  Falny-IVrot  pattern  showing  the  effort  of  our  reflection 
iiimli'  -elect or  ( RMS)  when  the  ruhy  is  tc-niiN-rntiirc  tuned  to 
resonance. 


I  in.  o.  I  :tht\  -I1'  t"i  i ii  -Iiouimu  tin  etiVci  hi  adding  ami  tun-  Fig.  0,  Fabry-I Yrot  pattern  showing  the  effect  of  adding  our  crv|>- 
ing  on i  iian-mi--ion  iiiim|<  -i  lector  IMS)  to  the  UMS.  tocyanine  dye  mode  selector  (CMS)  to  the  TMS  and  RMS. 


Fin  7.  Faliry-l’erot  pattern  showing  the  effect  of  adding  and  tun-  Fig.  8.  Fabry-I’erot  pattern  showing  the  effect  of  adding  our  CMS 
ing  our  TMS  to  the  RMS.  to  the  TMS  and  RMS. 
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width  for  the  ruby  laser  with  the  instrument  resolution 


indicates  thnt  the  corrected  ruby  linewidth  is  about  200 
Me.  This  linewidth  is  considerably  greater  than  the  ap¬ 
proximately  10-Mc  transform  from  the  pulse  width  meas¬ 
urements.  This  problem  has  not  been  reconciled,  and 
possibly  future  work  will  aid  in  our  understanding.  Some 
imssiblc  explanations  could  be  the  frequency  spreading 
associated  with  the  presence  of  a  large  number  of  trans¬ 
verse  modes,  each  having  a  slightly  different  frequency 
or  weak  components  separated  from  the  main  one  by 
the  200-Me-mode  interval  of  the  laser.  Also,  there  prob¬ 
ably  arc  independent  regions  with  slightly  different  path 
lengths  giving  rise  to  slight  frequency  spread.  Both  of 
these  hypotheses  are  consistent  with  our  initial  observa¬ 
tions  with  a  transverse-mode-controlled  laser.  Here  the 
line-width  of  individual  modes  appears  instrument- 
limited,  but  selected  portions  of  the  ruby  are  used  for 
laser  action. 
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MODE-STRUCTURE  INDEPENDENCE  OF  STIMULATFD 
RAMAN-SCATTERING  CONVERSION  EFFICIENCIES* 

F.  J.  McClung,  W.  G.  Wagner,  and  D.  Weinert 

Hughes  Research  Laboratories,  Malibu,  California 
(Received  13  April  1965) 


The  power  density  needed  to  convert  a  giv¬ 
en  amount  of  laser  energy  to  stimulated  Raman- 
scattered  radiation  has  been  found  to  be  about 
an  order  of  magnitude  less  in  nitrobenzene  than 
theoretically  expected  for  the  case  of  a  Raman 
cell  external  to  the  laser  cavity.1  It  has  been 
suggested1  ,z  that  the  multimode  character  of 
the  laser  pump  might  be  responsible  for  this 
disagreement.  Bloembergen  and  Shen2  have 
given  a  semiquantitative  estimate  of  the  enhance¬ 
ment  of  the  Raman  gain.  It  indicates  that  with 
a  typical  multimode  laser  a  Raman  gain  *-4  to 
8  times  greater  than  that  for  a  single-mode 
laser  should  be  produced.  It  is  the  purpose 
of  this  communication  to  show  that  the  Raman 
gain  in  nitrobenzene  for  ..  single-mtJe  laser 
pump  is  precisely  that  for  a  multimode  laser 
pump,  and  to  suggest  other  reasons  for  the 
above-mentioned  disagreement. 

In  this  experiment,  the  collimated  beam  from 
our  giant  pulse  laser  was  directed  onto  a  1.45- 
mm  aperture  70  cm  from  the  laser,  and  the 
portion  of  the  beam  transmitted  by  the  aperture 
passed  through  a  10- cm  cell  of  nitrobenzene 
onto  a  MgO  diffuse  reflector  placed  50  cm  af¬ 
ter  the  cell.  The  power  of  the  laser  and  the 
first  Stokes  line  was  monitored  by  suitably 
filtered  fast  photodetectors  that  sampled  the 
diffuse  reflection  from  the  MgO.  The  cell  was 
tilted  at  -3°  with  respect  to  the  laser  beam. 

This  arrangement  is  like  that  described  in 
more  detail  by  Weiner,  Schwarz,  and  McClung,1 
except  that  the  aperture  and  cell  are  now  fur¬ 
ther  from  the  laser  and  closer  to  ii*e  MgO. 

Also,  another  beam  splitter  to  allow  a  measure¬ 
ment  of  the  far-field  pattern  has  been  added. 

Our  laser  has  been  mode-selected  to  give 
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single-transverse  and  longitudinal  mode  be¬ 
havior.  The  details  of  the  mode -selection  tec*'  ■ 
niques  will  be  given  elsewhere.*  The  transvers* 
mode  structure  was  determined  with  the  aid 
of  a  1-m  focal-length  camera.  The  longitudi¬ 
nal  mode  structure  was  determined  with  the 
aid  of  a  2-cm  spaced  Fabry-Perot  etalon  with 
A/80  flat  plates  of  transmission.  When 
completely  mode-selected,  the  laser  produced 
2  MW  of  power  in  a  beam  whose  divergence 
equaled  the  diffraction  limit  corresponding  to 
the  laser-beam  diameter.  When  not  mode-se¬ 
lected,  the  laser  output  was  '10  MW,  the  beam 
divergence  ~1.5xi0"s  rad,  and  the  spectral 
width  cm”1.  The  pulse  length  was  -30  nsec 
for  both  cases.  The  results  of  the  conversion- 
efficiency  measurements  for  mode-selected 
and  non-mode-selected  lasers  are  shown  in 
Table  I.  The  experimental  arrangement  is  the 
same  for  both  cases.  The  relative  error  for 
the  power  measurements  is  -5%.  Measurements 
at  other  power  densities  for  our  mode-selected 
laser  gave  a  conversion-efficiency  curve  which 
agreed  very  well  with  our  previous  curve  for 
a  non-mode-selected  case.1  The  data  of  Ta¬ 
ble  I  indicate  strongly  that  the  anomalously 
high  gain  is  not  caused  by  the  laser  mode  struc¬ 
ture. 

The  theoretical  gain  in  nitrobenzene  at  20 
MW/cm  is  0.028  cm”1.  This  gain  is  computed 
using  a  formula  of  Hellwarth*  and  a  recently 
measured  peak  Raman-scattering  cross  sec¬ 
tion  of  1.3  i  0.4  cm"2.  This  cross  section  agrees 
within  the  expected  error  with  our  earlier* 
and  less  accurate  measurement  of  2.3  a  1.2 
cm"1.  It  also  agrees  with  the  recent  measure¬ 
ment  of  Damen,  Leite,  and  Porto*  in  the  fol- 
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Table  I.  Conversion  efficiency  for  mode-selected  and  non-mode-sei 

-cted  pumping  of  nitrobenzene. 

Description  of 
mode  structure 

Incident-peak 
power  density 
(MW/cm*) 

Power  converted  to 
first  Stokes 
radiation 
<%) 

Two  7 EM00  modes  of  -4:1 
intensity  ratio  separated 
by  ~^00  Mc/sec 

19 

0.6 

-100  transverse  modes 
-10-20  longitudinal  modes 

19.5 

0.6 

lowing  sense.  If  we  scale  our  data  for  benzene 
cross  sections  by  the  ratio  of  the  nitrobenzene 
cross  sections  (1.3/2. 3)  and  then  compute  the 
cross  section/molecule  at  6328  A,  we  find  a 
vaiue  of  0.7±0.2x  10“**  cm*.  Here  we  assume 
a  A*  wavelength  dependence  and  use  the  observed 
angular  dependence  for  Raman  scattering. 

This  value  agrees  (within  the  stated  error) 
with  their*  value  of  0.56*0.1  xiO“21  cm2. 

To  estimate  the  Raman  gain  at  20  MW/cm2 
from  our  conversion-efficiency  data,  we  as¬ 
sume  noise  is  amplified  during  a  double  pass 
(through  the  cell  to  the  laser  reflectors  and 
back  through  the  cell)  to  the  observed  level 
of  10*  W/cm*.  The  dominant  contribution  to 
the  noise  which  initiates  the  generation  of  the 
Raman-shifted  radiation  comes  from  the  quan¬ 
tum-mechanical  zero-point  vibrations  of  the 
electromagnetic  field.  In  a  spectral  interval 
of  dK,  and  within  a  cone  of  opening  angle  9, 
centered  on  the  axis  of  the  laser  beam,  there 
is  a  zero-point  power  flux7  given  by  (ir0*JW 
nx*)(dx/x).  Our  measurements  indicate  mat 
the  bulk  of  the  Stokes  radiation  emerges  with¬ 
in  an  angle,  9,  which  is  not  greater  than  0.017 
rad,  and  that  the  fractional  wavelength  spread, 
rfx/x,  is  of  order  3x10“*,  and  thus  the  start¬ 
ing  powe>  flux  is  estimated  to  be  "1.2  x  10“’ 
W/cm*.  Then  10*  W/cm*  =  1.2x  10”*  W/cm2 
xexp(ffx20  cm),  so  that  the  gain  per  cm  *  =  0.80 
cm”1.  The  gain  is  only  slightly  lower  if  cal¬ 
culated  for  multiple  double  passes  where  the 


feedbav.it  is  supplied  by  diffuse  reflection  from 
the  cell  window  and  the  MgO  reflector,  sinre 
the  solid  angle  of  acceptance  imposed  by  our 
geometry  is  so  small.  Thus  the  observed  gain 
is  at  least  25  times  greater  than  that  calculated 
from  the  cross  sections. 

We  are  very  grateful  for  the  penetrating, 
stimulating,  and  helpful  comments  of  Dr.  Shaul 
Yatsiv  during  the  course  of  this  work. 
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7It  is  assumed  that  the  Raman  radiation  is  strong: \ 
polarized,  in  accord  with  the  experiments.  The  wave¬ 
length  to  be  used  in  the  formula  is  that  of  the  Stokes 
radiation,  inside  the  Raman  cell,  and  n  is  the  linear 
index  of  refraction  at  the  Stokes  frequency  w. 
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The  Self-Focusing  of  Light  of  Different  Polarizations^ 

D.  H.  Close,  C.  R.  Giuliano,  R.  W.  Hellwarth, 

L.  D.  Hess,  F.  J.  McClung,  and  W.  G.  Wagner 

Hughes  Research  Laboratories,  Malibu,  California 

Abstract  —  In  an  effort  to  verify  the  hypothesis  that  molecular  reorien¬ 
tation  (ac  Kerr  effect)  is  responsible  for  the  self-focusing  of  light  beams 
in  certain  liquids,  we  have  measured  and  compared  the  thresholds  and 
other  characteristics  of  self-focusing  for  circularly  and  linearly  polarized 
beams  incident  on  these  liquids.  We  show  that  for  plane  waves  propagat¬ 
ing  in  a  homogeneous,  isotropic,  ensemble  of  molecules  having  anisotropic 
polarizability  tensors,  the  non-linear  index  should  be  four  times  as  great 
for  linearly  as  for  circularly  polarized  waves.  The  hope  that  this  differ¬ 
ence  in  indices  would  be  reflected  in  a  four  fold  increase  in  the  threshold 
pover  for  self -trapping  when  circularly  instead  of  linearly  polarized  light 
is  incident  was  not  realized.  In  practice,  the  increase  was  always  found 
to  be  much  less.  However,  in  every  case  studied,  the  trapped  light  from 
a  beam,  circularly  polarized  to  better  than  one  part  in  200,  was  markedly, 
if  not  completely,  depolarized  as  soon  as  self -trapping  could  be  detected. 

'  This  work  has  been  supported  in  part  as  a  part  of  Project  DEFENDER 
under  the  joint  sponsorship  of  the  Advanced  Research  Projects  Agency, 
the  Office  of  Naval  Research,  and  the  Department  of  Defense. 
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We  show  qualitatively  when  this  polarization  instability  should  exist  for 
all  but  linearly  polarized  light  and  for  a  variety  of  nonlinear  mecha¬ 
nisms.  However,  in  the  absence  of  even  an  approximate  quantitative 
theory  of  the  self -trapping  of  light  that  is  not  linearly  polarized,  the 
comparative  measurements  of  thresholds  cannot  be  said  either  to 
verify  or  disprove  the  hypothesis  of  molecular  reorientation. 

I.  INTRODUCTION 

Self-focuring  and  trapping  of  high  power  electromagnetic  beams 

1  2 

was  foreseen  by  Askar'yan,  Talanov,  and  Chiao,  Garmire,  and  Townes, 

4-8 

and  has  now'  been  observed  in  many  liquids.  In  this  effect  the  third 

order  nonlinear  index  of  refraction  of  a  substance  causes  the  incident 

light  beam  to  focus  or  "trap"  itself  into  one  or  more  narrow  beams  or 

"filaments.1  "  In  their  original  discussion  of  self-focusing  at  optical 
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frequencies,  Chiao,  et  al .  ,  estimated  the  relative  thresholds  for  self¬ 
trapping  that  might  be  expected  from  each  of  three  physical  mechanisms: 
molecular  reorientation  in  the  strong  optical  fields  (ac  Kerr  effect), 
electrostriction,  and  the  nonlinear  electronic  polarizability.  Recent 
comparisons  of  these  f‘°  imated  thresholds  with  thresholds  observed  in 
liquids  which  readily  exhibit  stimulated  Raman  scattering  have  suggested 
that  molecular  reorientation  is  most  likely  the  mechanism  underlying  the 

self-trapping  (a^d  hence  also  the  anomalously  low  Raman  and  Brillouin 
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thresholds)  in  these  liquids.  In  this  paper  we  report  the  results  of 
attempts  to  verify  this  molecular  reorientation  hypothesis  by  measuring 
the  ratio  of  the  thresholds  for  the  self-trapping  of  a  linearly  polarized 
beam  to  the  threshold  for  a  circularly  polarized  beam.  We  have  found 
an  instability  in  the  nonlinear  propagation  of  circularly  polarized  light 
which  has  made  interpretation  of  the  threshold  results  difficult. 

Both  a  linearly  polarized  plane  wave  and  a  circularly  polarized 
plane  wave  are  solutions  of  the  nonlinear  Maxwell  equations  for  a 
spatially  local  nonlinear  polarizability  proportional  to  the  field  cubed. 

In  general,  a  linearly  Dolarized  wave  propagates  with  a  different  (non¬ 
linear)  index  of  refraction  than  a  circularly  polarized  wave.  If  the 


m 
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nonlinear  index  arises  from  molecular  reorientation,  it  is  a  straight¬ 
forward  matter  to  show  that  the  real  nonlinear  part  of  the  index  for 
linear  polarization  should  be  four  times  that  for  circular,  regardless 
of  the  shape  of  the  molecules  being  oriented.  If  one  then  supposes 
that  at  least  the  onset  of  trapping  for  each  case  can  be  dealt  with  by 

geometrical  optics  incorporating  nonlinear  indices,  as  has  been  done 
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by  Kelley,  the  power  required  to  begin  trapping  a  circularly  polar¬ 
ized  beam  should  be  four  times  that  required  to  self -trap  a  linearly 
polarized  beam. 

Our  experiments  show  that  only  between  one  and  three  times 
as  much  circularly  polarized  power  is  required  in  the  many  liquids 
we  have  examined.  However,  we  have  found  that,  although  light 
trapped  from  a  linearly  polarized  input  beam  emerges  linearly 
polarized,  light  trapped  from  a  circularly  polarized  input  beam  is 
largely,  if  not  completely,  depolarized  as  soon  as  we  have  been  able 
to  observe  it.  This  fact  we  attribute  to  an  instability  in  the  propaga¬ 
tion  of  circularly  polarized  light.  We  show  that  this  instability  may 
exist  whenever  the  coefficient  of  the  term  (E  •  E  )E  in  the  expan - 

__  WWW 

sion  of  the  polarization  amplitude  P^  is  positive.  (E^  is  the  ampli¬ 
tude  of  the  electric  field  at  frequency  w.  )  This  coefficient  has  been 
called  B/2  by  Maker,  et^al.  ,  who  measured  it  in  a  variety  of  liquids 
and  found  it  always  to  be  positive.  The  coefficient  B  will  be  positive 
if  molecular  reorientation  or  electronic  nonlinearities  dominate  and 
doubtless  will  be  positive  for  many  other  possible  physical  mechanisms 
as  well.  We  have  not  discovered  any  mechanism  for  which  the  coeffi¬ 
cient  should  be  negative  (where  no  other  frequencies  are  present  in 
strength),  but  it  should  be  zero  for  stationary  electrostriction.  The 
fact  that  thresholds  for  trapping  do  not  agree  with  the  simple  predictions 
of  molecular  reorientation  does  not  rule  out  molecular  orientation  as  the 
main  mechanism  because  the  nonlinear  propagation  of  circularly  polarized 
light  appears  to  be  incorrectly  treated  by  all  the  techniques  devised  to  d'te. 
Only  the  nonlinear  index  for  linearly  polarized  light  can  be  obtained  from 
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trapping  measurements  in  liquids  studied  to  date,  for  a  linear  state 
of  polarization  is  the  only  state  which  propagates  in  a  manner  which 
can  be  approximated  by  the  scalar  nonlinear  equations  studied  so  far. 


II.  THEORY 


AC  Kerr  Effect 


In  the  present  experiments  we  are  interested  in  situations 
where  the  thermal  energy  kT  is  much  larger  than  a  rotational 
quantum  of  energy,  which  is  in  turn  much  larger  than  the  interac¬ 
tion  energy  of  a  molecule  with  the  strong  electric  fields.  We  are 
therefore  led  to  describe  the  molecules  in  the  liquids  classically. 
The  probability  P(J2)  for  a  molecule  to  have  its  major  dielectric 
axes  oriented  at  the  Euler  angles  0,  \\i  (symbolized  by  12)  is 

given  by 


P(«) 


exp  [-U(f2)/kT] 

7®  e^TOJA'T 


(1) 


where  the  energy  U(J2)  of  a  molecule  interacting  with  the  field  is 
given  by 


U(£l)  =  -  j  S._.  F.(t)  f.(t)  a^W)  ,  (2) 

the  bar  indicating  a  time  average  over  many  cycles  of  the  local  optical 
electric  field  F^(t).  The  indices  i,j  represent  the  spatial  indices 
x,  y,  z.  Here  a^(I2)  is  the  polarizability  tensor  for  a  molecule  oriented 
at  the  angles  J2.  When  the  principal  a::es  of  the  dielectric  ellipsoid  are 
oriented  along  the  '£,  y,  z  axe's  then  is  diagonal,  its  three  elements 
having  the  values  a,b,  and  c. 
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In  order  to  compute  the  desired  nonlinear  dielectric  constant, 
we  need  only  the  statistical  average  <CLj>  of  over  angles  ft  to 

second  order  in  the  fields  F^(t).  It  is  a  straightforward  matter  to 
perform  the  averaging  with  the  weight  function  (1)  and  obtain 

<a.  >  -  a  6. .  +  y. .  (3) 

iJ  ij 

where  a  =  (a  +  b  +  c)/3  is  the  usual  linear  optical  polarizability  of  a 
molecule  and 


(a 


b)2  +  (b  -  c)2  +  (a  -  c)2 

- 90TT - 


■ij  J  90  kT  |  2k,l  (36ik  6j I  '  6ij  6kf}  *k(t)Ff (t) 

(4) 

is  the  nonlinear  part  <  f  the  polarizability.  Although  for  a  single  plane 
wave  only  the  parts  of  y^  having  i,  j  =  x,  y  are  needed,  for  a  general 
treatment  of  trapping  all  space  components  of  y^  would  be  needed.  We 
may  take  the  local  field  F.(t)  to  be  related  to  the  macroscopic  electric 

renz  formula,  provided  that  the 
Then,  if  p  is  the  number  density 


field  E^(t)  by  the  linear  Lorentz- Lorenz  formula,  provided  that  the 
nonlinear  change  in  index  is  small 
of  molecules, 


F.(t)  =  E.(t)/(1  -  4irap/3)  .  (5) 


Following  further  application  of  the  Lorentz-Lorenz  local  field 
approximation,  the  polarizability  tensor  of  (4)  predicts  that,  for  a 
linearly  polarized  plane  wave  of  peak  amplitude  Eq,  the  linear  index 
of  refraction  n  is  altered  by  an  increment  An^  given  by 


An.  .  ((a-b,2t(b  -c)f  +(c-a)2\/^  wpEa 


(6) 


i  -  \  im? t 

For  a  circularly  polarized  plane  wave,  (4)  leads  to  an  increment  An 
to  the  linear  index  n  which  is  related  to  the  increment  in  (6)  by 


An  =  An .  /4  . 

c  I 


(7) 
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Many  other  relations  are  contained  in  the  expression  (4)  for 
the  nonlinear  polarizability  including  the  usual  expression  for  the 
ac  and  dc  Kerr  constants  (the  latter  for  nonpolar  molecules  only) 
which  follow  from  the  relation  of  Langevin  between  the  change  in 
index  An^  for  light  polarized  parallel  to  a  strong  linearly  polarized 
field  and  that  An^  for  light  polarized  perpendicular  to  it: 


An 


ll  =  "  Ani  /2 


(8) 


Stability  of  Polarization  in  Trapping 

If  an  electric  field  having  space  components  of  the  form 

Re  E.(£)e  la*  exists  in  an  isotropic  homogeneous  medium  with  which 

^  th 

it  interacts  locally  in  space,  then  the  i  spatial  component  of  the 

*■  icjt 

polarization  density  may  be  written  Re  P^(£)e"  where 

pi  ■  *jej  +  V,  W?  (9) 


to  third  order  in  the  fields.  From  the  isotropy  of  the  medium,  it  follows 
that  we  may  write 


Xijki  =A  Vw  +  7  B6il6jk  <10> 

where  we  have  followed  the  nomenclature  of  Maker,  et  al.  jn  a  jOS8_ 
less  medium  with  a  nonlinear  susceptibility  of  the  form  (10),  Maxwell's 
equations  are  satisfied  by  (among  many  other  things)  two  parallel  cir¬ 
cularly  polarized  plane  waves  propagating  with  indices  altered  by 
increments 


6  n±  =  (2  ir/n)  [a  E±E^  \  (A  +  B)  E^E*] 


(11) 
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respectively  for  the  right  and  left  circularly  polarized  waves  of  ampli¬ 
tudes  E±  =  (E^  ±  i  E^)/\T2?  However,  the  propagation  of  a  mono¬ 

chromatic,  finite  size,  diffraction  limited  parallel  beam  cannot  be 
simply  described  by  several  indices  as  in  (11).  Nevertheless,  we  can 
guess  what  might  happen  to  such  a  finite  beam,  launched  in  an  arbitrary 
state  of  polarization  from  a  plane  at  z  =  0,  by  a  study  of  (11).  For 
some  region  of  z  >  0  one  expects  that  the  circular  components  will 
essentially  just  have  their  phase  relations  altered  according  to  the 
difference  in  indices  predicted  by  (11).  -That  is-,  the  major  axis  of  the 
initial  polarization  ellipse  will  be  rotated  by  an  angle  p  tportional  to 
Bz.^  However,  beyond  a  certain  distance  zq,  the  fact  that  the  radial 
index  change  (caused  by  the  variation  in  intensity  across  the  beam)  is 
different  for  one  circular  component  than  for  the  other  (unless  |E+|  = 

|  E  |)  suggests  that  one  component  will  be  self-focused  more  strongly 
than  the  other,  and  hence  the  polarization  in  the  trapped  parts  of  the 
beam  will  change  with  length. 

If  A  and  B  are  positive  (which  seems  to  be  the  case  in  all 
materials  studied  so  far  and  which  holds  for  the  a*  Kerr  effect),  then 
the  weaker  of  the  two  circular  components  is  seen  to  be  self -trapped 
before  :he  stronger  because  of  its  greater  increase  in  index.  It  is  this 
tendency  of  the  weaker  cirulcar  component  to  increase  in  intensity  that 
we  believe  underlies  the  observed  depolarization  of  the  filaments  formed 
from  a  nearly  circularly  polarized  beam.  A  perfectly  circularly  polar¬ 
ized  incident  beam  is  impossible  to  obtain  and  hence  any  nominally  cir¬ 
cularly  polarized  beam  can  be  expected  to  degenerate  into  randomly 
oriented  linearly  polarized  filaments.  If  the  incident  beam  is  linearly 

^It  was  by  measuring  the  rotation  of  the  polarization  ellipse  that 
Maker,  et  al.  ,  originally  attempted  to  measure  the  B  coefficients 
for  various- materials.  10 


polarized,  then  any  slight  deviation  of  equality  between  E+  and  E 
will  tend  to  trap  the  weaker  component  and  restore  the  perfectly 
linear  polarization.  Hence,  only  linearly  polarized  light  can  be 
thought  of  as  a  single  stable  propagating  component  and  be  treated 
approximately  with  scalar  nonlinear  equations.  It  is  to  be  empha¬ 
sized  that  no  approximate  quantitative  treatment  of  the  self -focusing 
of  other  than  a  linearly  polarized  beam  has  been  put  forth  as  yet, 
and  the  arguments  presented  above  on  the  basis  of  (11)  are  qualitative 
at  best. 

We  note  from  (4)  that  one  may  easily  conclude  that  B/A  =  6 

for  the  ac  Kerr  effect.  For  a  nonlinear  electronic  polarizability  it 

is  easy  to  show  that  A  =  B;  and  for  static  electrostriction  that  B  =  0. 

A  weak  left  handed  component  in  the  presence  of  a  strong  right  handed 

one  is  focused  at  a  power  (2  A  +  B)/2  (A  +  3)  times  the  threshold  to 

3  t 

trap  linearly  polarized  beam.  ' 


The  birefringence  induced  inawe?.kbeam  of  light  at  5000  %.  by  a  strong 
linearly  polarized  beam  at  6943  %  has  been  measured  in  several  liquids 
by  Guy  Mayer  and  Francois  Gires  in  "Action  d'une  onde  lumineuse 
intense  sur  l'indice  de  retraction  des  liquids,  "  Comptes  Rendus,  Vol. 
258,  pp.  2039-2042,  February  1964.  The  coefficient  which  they  used 
to  describe  their  results  (and  called  B0  in  analogy  with  the  dc  Kerr 
coefficient),  is,  in  general,  independent  of  the  coefficients  A  and  B 
discussed  above,  as  B  describes  the  effect  of  two  frequencies  separ¬ 
ated  by  more  than  the  thermal  relaxation  rates.  However,  if  molecu¬ 
lar  reorientation  were  entirely  responsible  for  the  induced  birefringence 
and  self-trapping,  and  if  the  two  frequencies  in  the  former  are  not  too 
far  apart,  it  is  easy  to  deduce  from  our  eq.  (4)  that  the  B0  and  B  would 
be  related  by  B  =  n\B0/4ir  where  \  is  the  free  space  wavelength  of 
the  weak  birefringent  beam.  Using  this  relation  to  deduce  values  for 
B  from  the  measurements  of  BQ  by  Mayer  and  Gires  gives  values 
roughly  an  order  of  magnitude  higher  than  those  obtained  in  Ref.  10. 
However,  Wang  and  Racette  have  recently  remeasured  B  using  unfocused 
rather  than  fc'.used  beams  and  reported  (at  the  1966  International 
Quantum  Electronics  Conference)  much  higher  values  than  were  reported 
in  Ref.  10.  Using  these  new  B  values  and  some  refined  measurements 
of  Bq  reported  by  Mayer  at  the  same  conference,  we  calculate  the  ratio 
R  =  n\B0/(4trB)  to  be  1.  3,  1.  3,  1.  9,  and  1.  9,  respectively,  for  CS2» 
nitrobenzene,  toluene,  and  benzene.  These  remaining  deviations  of  R 
from  unity  cast  some  doubt  on  the  molecular  reorientation  hypothesis, 
but  may  still  be  the  result  of  experimental  uncertainties. 
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III.  EXPERIMENTAL  RESULTS 


We  have  measured  the  relative  thresholds  for  self -trapping  of 

linearly  and  circularly  polarized  light  in  a  number  of  liquids  by  two 

methods.  The  first  method  was  to  measure  photometrically  the 

thresholds  for  stimulated  Raman  scattering.  This  gives  a  fairly 

accurate  measure  for  the  trapping  thresholds  in  those  liquids  which 

show  strong  self-focusing  and  Raman  effect;  the  thresholds  for  the 
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appearance  of  both  then  coincide.  The  second  method  is  to  look 
photographically  for  the  first  appearance  of  self-focusing  in  the  time 
averaged  beam  profile  at  the  exit  face  of  the  liquid  under  study.  In 
both  types  of  experiments  linearly  polarized  light  was  found  to  produce 
like -polarized  filaments  and  Stokes  output,  whereas  circularly  polarized 
light  produced  depolarized  filaments  and  Stokes  output. 

Photometric  Measurements 

For  the  photometric  measurements,  the  unfocused  beam  of  a 
ruby  laser  was  first  passed  through  crossed  Gian  polarizers  which 
served  as  a  variable  attenuator  and  supplied  a  purely  linearly  polar¬ 
ized  output  beam.  The  unfocused  beam  then  passed  through  a  10  cm 
cell  containing  the  liquid  in  which  self -trapping  was  to  be  studied.  A 
portion  of  the  output  beam  from  this  cell  was  split  off  and  monitored  by 
silicon  photodiodes  for  Stokes  energy  content  after  passing  through  a 
Corning  5-56  filter,  two  Wratten  89B  filters,  and  appropriate  attenuators. 
When  the  liquid  cell  was  removed,  the  laser  produced  no  signal  at  the 
photodetector,  thus  assuring  that  only  Stokes  light  was  being  measured. 

Typical  results  for  the  Stokes  output  powers  plotted  as  a  function 
of  either  linearly  or  circularly  polarized  inputs  are  given  for  benzene, 
nitrobenzene,  and  CS2  in  Fig.  1.  From  such  data  the  trapping  threshold 
ratios  (linear  to  circular)  for  various  liquids  is  given  in  Table  I.  These 
ratios  are  seen  to  be  consistently  well  below  the  value  4  expected  if  cir¬ 
cularly  polarized  light  would  produce  circularly  polarized  filaments  via 
the  ac  Kerr  effect. 
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Fig.  I.  Plot-s  of  the  Raman  Stokes  output  energy  versus 
ruby  laser  energy  incident  on  a  10  cm  cell  of  (a) 
nitrobenzene,  (b)  carbon  disulphide,  and  (c) 
benzene.  Triangular  points  represent  linearly 
polarized  incident  light,  circular  points  represen' 
circularly  polarized  light.  The  zero  of  the  log 
Stokes  output  energy  corresponds  roughly  to  a 
conversion  of  10"^  of  the  laser  energy  to  Stoke 
energy.  The  incident  energy  was  not  calibrated 
absolutely.  The  dotted  lines  represent  pulses 
producing  no  detectable  Stokes  output. 


TABLE  I 


Ratios  of  Incident  Light  Powers  Required  to  Reach  Threshold 
for  Stimulated  Reiman  Scattering  (hence  self-focusing)  in 

Various  Liquids 


Liquid 

Threshold  Ratio 
Circular:  Linear 

Ni  robenzene 

1  ■  4(2) a 

Ca*  bon  Disulfide 

2.0(3) 

Benzene 

2.  2(?) 

B  romobenzene 

1-7(5) 

Toluene 

1.2(3) 

Benzoyl  Chloride 

1.3(4) 

Acetone 

1.3(5) 

The  figures  in  parentheses  give  the 
uncertainty  in  the  last  figure  estimated 
from  the  scatter  of  data  points. 

The  depolarization  of  the  Stokes  Raman  output  from  nitrobenzene 
was  also  studied  photometrically  by  imaging  the  exit  window  of  the  cell 
through  a  Fresnel  rhomb  and  a  Wollaston  prism  so  as  to  form  two  images 
on  a  boron  nitride  surface.  The  images  were  monitored  separately  by 
silicon  photodiodes  through  filters.  Depending  on  the  relative  orientation 
of  the  Wollaston  prism  to  the  Fresnel  rhomb,  signals  produced  by  the 
two  images  were  proportional  either  to  the  right  and  left  circular  compo¬ 
nent  intensities  or  to  the  intensities  of  the  components  linearly  polarized 
parallel  and  perpendicular  to  the  original  laser  polarization.  When  the 
cell  was  removed,  the  signals  showed  that  only  the  nominal  input  polariza¬ 
tion  (either  linear  or  circular)  was  present  to  within  one  part  in  200,  the 
limit  of  resolution.  When  Stokes  light  was  produced  by  linearly  polarized 
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input,  it  was  found  to  be  linearly  polarized  to  within  one  part  in  200; 
when  it  was  produced  by  a  circularly  polarized  input,  it  was  found  to 
contain  equal  intensities  of  both  circular  polarizations  and  of  both 
linear  polarizations.  These  results  applied  over  the  entire  power 
range  where  Stokes  light  could  be  observed  (rougly  seven  orders  of 
magnitude  variation  in  Stokes  output).  Also,  apertures,  from  1/2  to 
5  mm  diameter,  placed  in  the  laser  beam  did  not  affect  these  findings, 
even  when  the  cell  containing  the  liquid  was  placed  in  the  far  field  of  the 
1/2  mm  aperture. 

Photographic  Techniques 

In  the  photographic  measurements,  the  output  from  a  ruby  laser 
was  passed  through  attenuating  filters,  a  1.  5  mm  diameter  aperture  and 
a  rotatable  quarter  wave  plate  into  a  10  cm  cell  containing  the  liquid 
under  study.  The  output  end  of  the  cell  was  imaged  through  a  Fresnel 
rhomb  and  Wollaston  prism  onto  the  entrance  aperture  of  a  low -dispersion 
grating  spectrograph.  The  output  of  the  spectrograph  was  recorded  on 
Eastman  Kodak  IN  spectroscopic  plates.  A  single  laser  pulse  then  pro¬ 
duced  six  (or  eight)  separate  images  magnified  about  eight  times  and 
representing  two  orthogonal  polarization  components  of  the  incident  laser 
beam  and  the  first  and  second  (and,  in  the  case  of  CS^,  the  third)  Stokes 
shifted  output  beams.  The  Stokes  orders  all  appeared  ir;  filaments  which 
were  of  the  order  of  5  to  30  p  in  diameter,  several  times  smaller  than 
the  corresponding  filaments  (~30  -  50  p  diameter)  in  the  laser  beam.  The 
first  Stokes  filaments  were  of  the  same  order  of  intensity  as  the  laser 
beam  filaments. 

Ifthe  input  beam  was  linearly  polarized,  the  filaments  and  Stokes 

outputs  from  nitrobenzene  and  CS7  were  like -polarized  to  within  one 
3  C 

part  in  10  ,  the  limit  of  instrumental  resolution.  However,  if  the  input 
beam  was  circularly  polarized,  then  the  filaments  and  Stokes  outputs 
from  nitrobenzene  contained  equal  amounts  of  both  circular  polariza¬ 
tions  over  the  entire  range  of  observability  (~  three  orders  of  magnitude 
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in  the  first  Stokes  intensity).  The  right  and  left  circularly  polarized 
images  at  any  one  output  frequency  were  identical  in  every  detail. 

When  the  Wollaston  prism  was  rotated  to  test  for  orthogonal  linearly 
polarized  components  under  the  same  (circular)  input  conditions, 
the  total  intensity  of  one  polarization  in  all  of  the  filaments  at  a  given 
frequency  was  equal  to  the  total  filament  intensity  of  the  orthogonal 
polarization,  in  agreement  with  the  results  of  the  photometric  studies. 
However,  the  filament  pattern  for  one  linear  component  and  frequency 
did  not  now  agree  in  detail  with  the  pattern  for  the  orthogonal  polariza¬ 
tion  (at  the  same  frequency)  but  suggested  that  each  filament  might 
consist  mainly  of  light  linearly  polarized  in  some  random  direction. 

The  photographic  comparison  of  circular  and  linear  thresholds 
was  more  difficult  and  less  accurate  than  the  photometric  comparisons. 

A  plot  of  the  number  of  Stokes  filaments  observed  as  a  function  of  both 
linearly  and  circularly  polarized  input  energies  indicated  that  the  threshold 
ratios  were  less  than  four,  but  perhaps  somewhat  higher  than  those  in 
Table  I  for  CS2  and  nitrobenzene. 

Similar  experiments  with  circularly  polarized  light  trapped  in 
CS2  showed  some  slightly  different  results.  There  was  incomplete 
depolarization,  only  of  order  of  one  third  of  the  filament  energy  appear¬ 
ing  in  the  oppositely  polarized  laser  and  first  Stokes  filaments.  However, 
the  partition  of  energy  was  more  nearly  equal  in  the  second  Stokes  output 
and  became  essentially  equal  in  the  third  Stokes  output. 

In  the  abcve  and  previous  photographic  studies,^  in  which  the 
output  light  is  averaged  over  the  entire  (~25  nsec)  duration  of  the  laser 
pulse,  trapped  filaments  appeared  to  be  at  most  twice  as  intense  as  the 
average  beam  intensity,  and  hence  the  explanation  of  the  observed  Raman 
gains  (several  orders  of  magnitude  larger  than  predicted)  was  not  com¬ 
plete.  To  clarify  further  the  role  of  trapping  in  producing  anomalous 
Raman  gains,  we  made  time -resolved  photographs  looking  through  the 
spectrograph  at  the  exit  end  of  the  cell  (but  with  the  Fresnel  rhomb  and 
Wollaston  prism  removed).  These  showed,  as  expected,  that  at  their 
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peak  intensity,  trapped  fialments  were  one  to  two  orders  of  magnitude 
brighter  than  the  background.  The  laser  filaments  appeared  at  differ¬ 
ent  times  lasting  typically  3-6  nsec.  TI  2  first  Stokes  fialments  were 
observed  to  persist  for  a  shorter  duration,  ~1  to  2  nsec.  These 
observations  showed  qualitative  agreement  between  the  expected 
Rarran  gain  and  the  laser  intensity. 
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APPENDIX  IV 


Details  of  a  Cal  relation  of  the  Term  U22  in  the 
Fourth  Order  Electric  Free  Energy 

2  2 

Here  we  outline  the  calculation  of  =  <  <  u~  »  -  <  <  u~  »  . 

*  ?  ^33 

We  will  use  the  form  of  eq.  (37)  for  the  probability  N~L  p(£j  -  ^2)  ^  ^2^  ~2 

that  if  a  certain  molecule  is  in  d^r,,,  another  specified  molecule  is  in 

3  .  L 

d  ^2  •  We  will  use  Kirkwood's  superposition  approximation  for 

P3(Xr  X2'  ^3)  =  p3(  123)  [and  for  V‘4  p4(x1,  £ 2  ,  x3*  £4)  =  p4(1234)] 


-3 


which  is  the  probability  (per  unit  volume  for  each  molecule)  for  finding 
specified  molecules  at  £j,  £ ^ ,  £^  (and  £4)  simultaneously: 


p3(123)  ~  p(  12)  p(23)  p(  13) 


(Al) 


P4(  1234)  =  p(  12)  p(  1 3)  p(14>  p(23)  p(24)  p(34) 


(A2) 


Here  p(i2)  is  an  abbreviation  for  p(£19),  and  we  have  ignored  terms  of 
-2 

order  N  smaller  than  the  lead:ng  terms  of  (Al)  and  (A2).  (Terms  of 
order  N  ^  must  be  kept,  as  the  terms  in  U22  which  are  proportional 
to  will  cancel.  ) 

2 

The  integrals  in  4  <<  U2  »  fall  naturally  into  four  groups,  which 
we  call  Ij.I2.I3i  and  I4>  depending  on  whether  1  - ,  2 - ,  3 - ,  or  4-particle 
correlations  are  required: 

4  «*\  »  =  if14Ii  (A3) 

where 

II  =  z,  <<  (S'  ■  S?  ■  ■  a  •  S’)2  »  =  «  (A4) 


Equation  numbers  in  text  refer  to  those  in  the  main  body  of  the  report. 


A-  1 


h  ■-  2V «  (S' '  nr  ■  B 


■  ■  S  •  aV  •  E')2 


r^J  r*-> 


)  » 


(A5) 


V^pvy^'  •  ^  ■  a"’  S’ 


S'  ■  s?  ■  fiVY-  &Y  •  S'  » 


(A6) 


and 


V  Vpy^®'  '  a>1-  ^  ■  &v'  ®' 


S’  ■  /  •  BPv  ■  ft'  •  S’  » 


(A7) 


Here  the  summed  indices  are  always  unequal  in  any  I..  The  Ij  is  zero 
of  course  because  is  defined  to  be  zero.  The  angular  averages  in 

each  of  the  remaining  integrals  are  easily  performed  and  give 

I2  =  2a4E'4£Hiv[(l  +  A2/4)2  <(D^)2>  +  (A2/2  +  11  A4/l6)<Tr  '  & 

(A8) 


I.  =  4a4  E'4(l  +  A2)  £  <D,lvDvv> 
3  pv  y  xx  xx 


(A9) 


1.  =  a4  E,4£  a  <D,UV  D0V> 
4  pvpy  xx  xx 


(A10) 


Here  "Tr"  denotes  the  "trace"  or  sum  over  diagonal  space  elements,  and 
the  single  bracket  O  denotes  the  average  over  spatial  configurations. 

Using  the  form  (37)  for  the  two  particle  distribution  function  and 
(20)  for  D^V  gives  (neglecting  terms  ~l/N) 

tt  oo 

£  <(D>av)2>  =  N2V 

pv  xx' 


^  4tt  J  sin  v|j  dvj;  J  r2  dr(l  -  3  cos2i|i)2r  0 


=  Np  64  t r  / (45  t)  . 


A-  2 


(All) 


Similarly, 


TT 


CD 


^v< TrJJ^  ■  P^V> 


N2V  ^  1  Ztt  J  sin4)  J  r2  dr 


(1+3  cos^)r  ^  =  Np  32  it2/(3t)  . 


(A12) 


Therefore, 


I2  =  Nq4E'4  p2  4  tt2(32  +  136  A2  +  167  A4)/(45  Pt)  .  (A13) 


To  perform  the  average  in  (A9)  to  order  N,  we  write 


p(r)  =  (1  -  h(r))/<l  -  t/V) 


(A14) 


where  h(r)  ~  1  for  r  <  d  and  is  zero  otherwise.  Then,  with  (Al),  we 
obtain  to  order  N 


J  <D^V  Dvy>  =  N3V"2/d3  d3  15  D 
(ivy  xxxx  J  ^  2  ^3  12 


23  (1  '  h13J  ’ 


(Al  5) 


.12 


where  a  (1  -  h12),  hJ2  denotes  h(£2  -  £3).  etc.  The  integrals 

here  are  easily  performed  using  the  theorem  that  the  integral  J  of 
©12  over  any  volume  Q  is  related  to  an  integral  over  the  surface  S 
bounding  <u  by 


=  /d3*2D  12*-  / 


d2  £2  n  •  2  x12  •  2  rj2  ,  (Al 6) 


provided,  of  course,  that  r12  is  greater  than  the  molecular  diameter  d 
everywhere  on  the  surface  S.  (If  S  is  ellipsoidal  and  is  inside  it, 
recall  that  J  =  i.  If  S  is  spherical,  r^  is  outsiae  it  and  r^2  >  d 


A-3 


on  S,  then  J  =  QD^>  where  is  the  center  of  the  sphere  and  Q 

is  the  numerical  value  of  the  volume  bounded  by  S.  )  The  part  of  (A15) 

2  2 

independent  of  h.  ,  is  obviously  Np  t  .  With  the  aid  of  (A16),  the 

I J  2 

part  of  (A15)  involving  h^  integrates  to  Np  2tt  /3  from  which 

I3  =  4N  a4  E'4  p2(l  +  A2)  (I2  +  2tt2/3)  .  (A17) 

To  evaluate  1^  we  use  (A2)  to  obtain 

L  .  <D^V  D^V>  =  N(N  -  1)  (N  -  2)  (N  -  3)V"3  (1  -  t/V)'6 
pivp-y  xx  xx  '  '  '  ' 

I d3  *2  d3  X4  ®12  B34  (1  •  h13>  -  h14>  (1  ‘  h23>  ‘  h24> 

=  St=0Ki  •  (A18> 

The  right  hand  side  separates  naturally  into  terms  which  we  call 
where  i  =  0,  1,  ...  4  is  the  number  of  explicit  h -factors  in  the  inte¬ 
grand  of  the  term.  The  K  term  is  integrated  trivially  and  has  a  part 

2  °  2  2 

of  order  N  which  cancels  the  N  part  of  <<on7»  ,  and  a  part  of  order 

2  ^ 

N,  which  combines  with  «u,2»  from  (41)  (corrected  by  the  now 
important  factor  (1  -  t/V)  “ *)  to  give 

K  -  4  «u_»2  =  -  4  Np2  1 2(1  -  Pt)  .  (A19) 

O  u 

There  are  four  terms  in  K.  (there  being  four  possible  terms  linear  in 

1  2  3 

one  of  the  h's)  each  of  which  integrates  trivially  to  -f  Tp  giving 

K1  =  -  41  2p3  tN  .  (A20) 


A-4 


There  are  six  ways  to  form  integrands  quadratic  in  the  h's  in  (A18), 
four  of  which  have  one  space  point  as  an  argument  of  both  h's  and 
are  equivalent  to  each  other.  The  remaining  two  have  no  space  point 
as  arguments  of  both  h  factors  and  are  equivalent  to  each  other. 
Each  of  the  first  four  terms  is  seen  to  be  identically  zero  with  the 
aid  of  (A16),  and  the  last  two  can  be  straightforwardly  integrated 
with  (A  16)  to  give 


K2  =  0.  5066  tt2p3 tN  .  (A21) 

Almost  the  same  straightforward  integrations  are  required  for  each  of 
the  four  equivalent  terms  which  are  cubic  in  the  h's  and  one  finds  that 

K3  =  -  0.  11713  tt2p3tN  .  (A22) 

We  have  not  been  able  to  integrate  the  remaining  term  analyti 

cally.  We  estimate  that  it  is  an  order  of  magnitude  smaller  than 
because  the  four  h  factors  combine  to  limit  the  volume  accessible  to  the 
points  of  integration  much  more  than  in  K^.  Since  is  already  a 

minor  contribution  to  the  final  result,  we  feel  well  within  the  other  limits 
of  our  calculation  to  neglect  in  the  final  result: 

I4  -  4  «ii2»2  *  Np2( -41  2  +  3.  84  p t)  (A23) 

combining  all  the  above  terms  gives  the  U7~  of  (38). 
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